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1 Introduction 



1.1 Quark Mixing and CP Violation in the Standard Model 

In the Standard Model (SM), quark masses and their mixing and CP violation have a common 
origin. Local gauge invariance forbids bare mass terms for fermions, just as it forbids gauge 
boson masses. Masses are generated after spontaneous electroweak symmetry breaking, owing 
to the Yukawa couplings of the fermions to the Higgs condensate. Therefore CP violation and 
flavor-changing interactions also probe the electroweak scale, complementary to other studies at 
the energy frontier. Because the nature of the Higgs condensate is not known, it is interesting 
to study as many different aspects of it as possible. 

Historically, flavor physics and CP violation were excellent probes of new physics: the absence 
of the decay — > ^^^^ predicted the charm quark; CP violation in K^K^ mixing predicted 
the third generation; the K^K^ mass difference predicted the charm mass; the mass difference of 
the heavy and light mass eigenstates of the B^B^ mesons predicted the heavy top mass. These 
measurements put strong constraints on extensions of the SM, and imply that if there is new 
physics at the TeV scale, it must have a very special flavor and CP structure to satisfy these 
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constraints. 

The SM contains three sources of CP violation (if the quark masses are nonzero and nonde- 
generate). One of them occurs in the Cabibbo-Kobayashi-Maskawa (CKM) matrix |l|, ^ that 
describes the mixing of the quark generations, and the corresponding CP-violating parameter 
is of order unity. This source of CP violation is discussed in detail in this review. With the SM 
amended to include neutrino masses, CP violation is also possible in the mixing of leptons. We 
have no experimental information on the size of this effect yet. The third source of CP viola- 
tion occurs in flavor-conserving strong-interaction processes. However, the experimental upper 
bound on the neutron electric dipole moment implies that this {9QC]:,/16'ir'^)F^^Ff^'^ term in the 
SM Lagrangian is at best tiny, ^qcd ^ 10^^'^ (for a review, see Q). This is usually referred to as 
the strong CP problem. Thus we know that neither CP violation in the lepton sector nor that 
corresponding to ^qcd can play any role in quark-flavor physics experiments. 

Interestingly, there already exists strong evidence for CP violation beyond the SM. Assuming 
that the evolution of the universe began in a matter-antimatter symmetric state or that inflation 
washed out any previously existing asymmetry, the three sources of CP violation in the SM 
mentioned above cannot explain the observed baryon-to-entropy ratio, ub/s ~ 9 x 10^^^. The 
SM contains all three Sakharov conditions [Q] — baryon number violation, C and CP violation, 
and deviation from thermal equilibrium — required to generate such an asymmetry. However, 
the asymmetry generated by the SM is many orders of magnitudes too small [||. This implies 
that there must be CP violation beyond the SM, but it does not have to be in the quark 
sector, nor in flavor-changing processes.^ Nevertheless, because most new physics scenarios 
predict or naturally have observable effects in quark-flavor physics experiments, it is interesting 
to thoroughly test this sector of the SM. 



1.2 The CKM Matrix 

The Yukawa interaction of the quarks is given by 



Cy = -Yt, Q'u <t> d^R, - Yi; Qi, e 4>*u'^^ + h.c, (1) 

where y"''^ are 3x3 complex matrices, is the Higgs field, i and j are generation labels, and 
e is the 2x2 antisymmetric tensor. The left-handed quark doublets, and and 

are right-handed down- and up-type quark singlets, respectively, in the weak-eigenstate basis. 
When acquires a vacuum expectation value, (i?i>) = (0, f/\/2), Eq- (0) yields Dirac mass terms 
for quarks with 3x3 mass matrices 



^The baryon asymmetry may be explained by leptogenesis The discovery of neutrino masses makes the 
seesaw mechanism and the existence of very massive right-handed neutrinos likely. These may decay out of 
equilibrium through CP-violating interaction, generating nonzero lepton number, which can subsequently be 
converted by part into a baryon asymmetry by [B + L)-violating but [B — L)-conserving sphaleron processes 
present in the SM. 
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To move from the basis of the flavor eigenstates to the basis of the mass eigenstates, one performs 
the transformation 

jjuid)^uid)^u(d)^ = diag(m,(,),m,(,),mi(,)) , (3) 

where U^''^ and U^'^ are unitary matrices and the masses niq are real. The quark mass matrices 
are diagonalized by different transformations for the left-handed up- and down-quarks, which 
are part of the same SU(2)l doublet, 

By convention, we pulled out {U^)ij, so that the "misalignment" between the two transfor- 
mations operates on the down-type quark mass eigenstates. Thus the charged-current weak 
interaction is modified by the product of the diagonalizing matrices of the up- and down-type 
quark mass matrices, the so-called CKM matrix 

/ Ks Vub \ 

V = UlUf} = Vcd Vcs Vch ■ (5) 
V Vu Vts Vtb I 

However, the neutral-current part of the Lagrangian in the mass eigenstate basis remains un- 
changed, i.e., there are no flavor-changing neutral currents (FCNC) at tree level. 

Being the product of unitary matrices, V itself is unitary, VV'^ = I. This requirement and 
the freedom to arbitrarily choose the global phases of the quark fields reduce the initial nine 
unknown complex elements of V to three real numbers and one phase, where the latter accounts 
for CP violation. Because these four numbers effectively govern the rates of all tree- and loop- 
level electroweak transitions that involve the charged current, it is a compelling exercise to 
overconstrain V . If inconsistencies among different measurements occurred, it would reveal the 
existence of physics beyond the SM. 

The charged-current couplings among left-handed quark fields are proportional to the elements 
of V . The right-handed quarks have no Vl^-boson interaction in the SM, and the photon and 
gluon couplings are flavor-diagonal. For left-handed leptons, flavor-mixing can be introduced in 
a manner similar to that of quarks. 

The first CKM element, the Cabibbo quark-mixing angle 9c, was introduced in 1963 |j^ to 
explain the small weak-interaction decay rates for particles carrying strangeness. When CP 
violation was discovered in 1964 by the observation of the CP-odd decay K'^ 7r"*"7r~ |^, 
researchers had not yet perceived the intimate relation between the dynamical rules of quark- 
flavor mixing and the phenomenon of CP violation. Hence the terrain was open for speculations. 
In 1970, Glashow, Iliopoulos &: Maiani (GIM) Q used the unitary quark-mixing ansatz to 
postulate the existence of a fourth quark with quantum number charm to explain the observed 
suppression of strangeness-changing neutral currents (e.g., — > fi~^fi~). This mechanism yields 
the absence of tree-level FCNCs in the SM. In 1973 the concept of quark-flavor mixing and CP 
violation were unified when Kobayashi & Maskawa showed that for at least three generations of 
quarks, there would be enough physical degrees of freedom left in the quark-flavor mixing matrix 
to allow for a nonzero phase O]. This example demonstrates how, throughout the history of 
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particle physics, discoveries and developments in flavor physics have led to spectacular progress. 
The subsequent discovery of bottom and top quarks, and even a third lepton generation, as 
well as the observation of direct CP violation in the kaon system backed the KM idea. The B- 
factories provided the proof that the KM mechanism is the dominant source of CP violation at 
the electroweak scale with the measurement of sin2/3 in agreement with the KM expectation 

Beyond the description of CP violation, the three-generation CKM matrix provides a consis- 
tent framework for understanding a huge variety of observations. This includes understanding 
the very different neutral-meson mixing frequencies governed by FCNC box diagrams (the GIM 
mechanism is not effective for top quark loops and hence large mixing amplitudes are possible), 
to the precise predictions of rare decays of K and B mesons. Using the measurements of the 
K^K^ and B^B^ mixing frequencies, the KM ansatz could be used to predict the scales of charm 
and top-quark masses, respectively. Consequently, the hierarchical structure of the CKM matrix 
guides experimentalists in their quest for observables that can precisely determine its elements: 
a) The quarks occurring in the box diagrams govern the mixing frequency, and hence determine 
whether or not the neutral meson mixing can be observed experimentally; b) the lifetimes of 
the lightest mesons with a given flavor depend on the CKM elements in the dominant weak de- 
cays; to measure time-dependent CP asymmetries, the lifetime and the mixing frequency should 
be comparable; c) the CKM mixing hierarchy between the generations determines the relative 
size of the observable CP-violating effects, and which mixing and decay amplitudes contribute 
to the process. The conjunction of these considerations have led to the understanding that 
asymmetric-energy i3-factories would allow theoretically clean measurements of CP-violating 
observables in the Bd sector. As a result of the successful operation of the i?-factories, for 
the first time, significantly overconstraining the CKM matrix with many dynamically different 
processes, and hence testing the KM theory, was possible. This — and not the measurement of 
the CKM phase's value per se — has been and is the primary motivation for investing important 
resources into the physics of quark flavors. The overwhelming success of the KM theory hides 
however the problem that the origin of the observed flavor structure (mass hierarchy, mixing 
and CP violation) is not understood within the SM. 

1.3 Unitarity Constraints 

We have seen in the previous section that quark-flavor (or weak interaction) eigenstates and 
mass eigenstates are not equivalent. Hence the quark flavors mix through their couplings to 
the charged weak current. The CKM matrix, which describes the three-generation quark- flavor 
mixing in the SM, would be greatly simplified if, for instance, two of the up- or down-type 
quarks had equal masses and were therefore indistinguishable (see Sec. 1.3.2| ). Because there is 



no empirical evidence for such degeneracy, the full variable space for the CKM elements must be 
considered. Fortunately, the number of free parameters can be greatly reduced by very general 
considerations. Unitarity and the freedom to arbitrarily choose the global phase of a quark 
field, reduce the original 2n'^ unknowns (where = 3 is the number of generations) to (rig — 1)^ 
unknowns. Among these ng{ng — l)/2 are rotation angles and {ug — l){ng — 2) /2 phases describe 
CP violation. Three generations allow for only a single CP-violating phase. (In the lepton 
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sector and if neutrinos are Major ana fermions, there are Ug — 1 additional CP- violating phases, 
which are however not observable in oscillation experiments.) 

Among the infinite number of possibilities to parameterize the CKM matrix in terms of four 
independent parameters, we recall the two most popular ones for the purpose of this review. 

1.3.1 CKM Parameterizations 

Chau and Keung [|lO| proposed the "standard parameterization" of V. It is obtained by the 
product of three (complex) rotation matrices, where the rotations are characterized by the 
Euler angles 9i2, ^13 and 623, which are the mixing angles between the generations, and one 
overall phase 6, 

/ C12C13 S12C13 si3e~"^\ 

V = -S12C23 - Ci2S23Sl3e*'^ C12C23 - Sl2S23Sl3e*'^ S23C13 , (6) 
\ S12S23 - Cl2C23Sl3e*'^ -C12S23 - Sl2C23Sl3e*'' C23C13 / 

where Cij = cosOij and Sij = sinOij for i < j = 1, 2, 3. This parameterization satisfies exactly the 
unitarity relation.^ 

Following the observation of a hierarchy between the mixing angles, ,313 <^ S23 ■5i2 
Wolfenstein proposed an expansion of the CKM matrix in terms of the four parameters 
A, A, p and 77 (A ~ ~ 0.23 being the expansion parameter), which is widely used in the 
contemporary literature. We use the definitions to all orders 

S12 = A, 

S23 = , (7) 
sne-'^ = A\\p - iij) . 

Inserting the above definitions into Eq. (|6|) gives exact expressions for all CKM elements. Al- 
though expanding Eq. (|6|) in A is illustrative, 

/ 1 - ^A2 - iA^ A AX^{p-ir,)\ 

V=l -A l-iA2-iA4(l + 4^2) ^^2 +0(^5)^ 

\AX^{l-p- ir]) -A\^ + i^A^ (1 - 2{p + ii])) 1 - ^A^X'^ ) 

it is not necessary, nor is it necessary to truncate the expansion at any order in A. In Ref. 
and in this review no such approximation is made. 

1.3.2 The Jarlskog Invariant 

A phase-convention-independent measure of CP violation in the SM is given by 

Imdet ([M"M"^ M'^'M'^t]^ = 2 J {m^ - ml){ml - ml){ml - ml) 

X {ml - ml){ml - ml){ml - ml) . (9) 



^This phase 5 describes CP violation. It should not be confused with the CP-conserving hadronic phases that 
are introduced later. 
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Im {Vi.VklViVC^) = J E ^ikm Sjln , (10) 
m,n=l 

where Vij are the CKM matrix elements and £ikm is the total antisymmetric tensor. Hence 
owing to the unitarity of V, the nonvanishing imaginary parts of all quadri-products of CKM 
elements are equal up to their sign. One representation of Eq. (|^) reads, for instance, J = 
lm{VudVcsyus^*d) ■ ^ nonvanishing KM phase and hence CP violation requires J 7^ 0. The 
Jarlskog parameter expressed in the standard parameterization reads 

J = Cl2C23Ci3Si2S23Sl3Sin(5 , (11) 

and, using the Wolfenstein parameterization, 

J = A^X^T] (l - AV2) + 0(Ai°) « 3 X 10-5 . (12) 



The empirical value of J is small compared with its mathematical maximum of 1/(6^3) « 0.1, 
showing that CP violation is suppressed as a consequence of the strong hierarchy exhibited by 
the CKM matrix elements. Remarkably, to account for CP violation (see Eq. (^)) requires not 
only a nonzero J, but also nondegenerate masses of the up-type and down-type quarks. 

Physically meaningful quantities are phase-convention invariant. Such invariants are the 
moduli, \Vij\, and the quadri-products VijVkiy^u^kj (^-S-) Jarlskog invariant J). Higher- 
order invariants can be rewritten as functions of these [^. The Wolfenstein parameters A = 
I I / vTKjdP~+~n47F ^ = I^cft/Kisl/''^ are phase-convention invariant, whereas p and r] with 
p + ir] = V*i^/{AX^) are not. We use phase-invariant representations and formulas throughout 
this review. 

1.3.3 Unitarity Triangles 

The allowed region for the CKM phase can be elegantly displayed by means of the Unitarity 
Triangle (UT) described by the rescaled unitarity relation between the first column of the CKM 
matrix and the complex conjugate of the third column (i.e., corresponding to the meson 
system^) 

VM VcdV:, v^dV^ ^ ' ^ ^' ^ ' 

Note that twice the area of the non-rescaled UT equals the Jarlskog parameter J (this is true 
for any one of the six UTs, because they all have the same area). This provides a geometrical 
interpretation of the phase-convention invariance of J: A phase redefinition of the CKM matrix 
rotates the UT while leaving its area, and hence J, invariant. It is a remarkable property of 
the UT in the B system that its three sides are proportional to the same power of A (so all 



sides of the rescaled UT ( |13D are of order one), which predicts large observable CP violation. In 



comparison, the corresponding UT for the kaon sector is strongly flattened 



^The subscript in Bd is omitted in the following. 
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implying small observable CP violation. 

The UT (13) is sketched in Fig. || in the complex plane, where the apex is given by the 
following phase-convention independent definition fl^ , 



p + ir] 



(15) 



When the above equation is expressed in the Wolfenstein parameterization, one finds to all 
orders in A [Ol^ 



p + ir] 



VI - ^2A4(^ + vq) 



(16) 



The sides Ru and Rt of the UT (the third side along the real axis is normalized to unity) read 
to all orders, 



Ru 



VudV, 



ub 



p2 + rj'^ 



Rt 



The three angles of the UT, a, P, 7, are defined by 



VuV, 



tb 



1-pf+t- 



a 



arg 



VtdV, 



tb 



Vudv:, 



13 = arg 



cb 



vtdv:. 



7 = arg 



VudV,, 



ub 



Vcdv:, 



(17) 



(18) 



and the KM phase in the standard parameterization (^) is 5 = 7 + A^X^r] + ©(A^). 

At hadron machines there is copious production of mesons. Here the corresponding Uni- 
tarity Triangle, UT^, is obtained by multiplying the second column of the CKM matrix by 
the complex conjugate of the third column, which is equivalent to replacing all d-quarks with 
s-quarks. Appropriately rescaled, one has 







0(A2) + 1 + 0(1) 



The apex of the UT^ is defined accordingly by 



(19) 



(20) 



^Expanding Eq. ( |l5| ) in A leads to the well-known approximation 

p = p(l - aV2) + 0(A*) , 77 = ^7(1 - aV2) + OiX") . 
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which is a short vector of order A^. The smaU angle (3s opposing the short side is given by 

The quantity sin2/?s can be determined from a time-dependent analysis of Bg —>■ J/ip(l) decays. 



Af = {f\H\B) = J2Ake'^''e''^'=, Aj = {f\H\B) = J2Ake'^'= e-'^' . (22) 



1.4 CP Violation and Neutral-Meson Mixing 

The most precise and among the theoretically cleanest information about the phase of the 
CKM matrix stems at present from the measurements of time-dependent CP asymmetries in 
B decays, so we need to review briefly the relevant formalism. We use the B^ system as an 
example. In general, all forms of CP violation are related to interference phenomena, because CP 
violation is due to irreducible phases in the Lagrangian, which are observable only in interference 
experiments. 

The conceptually simplest form of CP violation, which can occur in both charged- and neutral- 
meson as well as baryon decays, is CP violation in decay. If at least two amplitudes with nonzero 
relative weak (cp^) and strong (6k) phases (that are odd and even under CP, respectively) 
contribute to a decay, 

k k 

Then it is possible that lAj/Ajl ^ 1, and thus CP symmetry is violated. If there are two 
contributing amplitudes then \Aj\'^ — \Af\'^ oc sin(0i — 02)sin(5i — ^2). Because these kinds 
of CP asymmetries depend on strong phases, their interpretation is usually model dependent 
(unless the strong phase can be eliminated by relating several decays to each other). The best 
established observations of this type of CP violation are Reej^ and the charge asymmetry in 

In neutral-meson decays there are other possibilities in which CP violation can occur. The 
two i?-meson mass eigenstates are linear combinations of the flavor eigenstates, 

\BL,H)=p\B'')±q\W), (23) 

with two complex parameters p and q. The CP symmetry is violated if the mass eigenstates are 
not equal to the CP eigenstates {CP violation in mixing). This happens if \q/p\ 7^ 1, i.e., if the 
physical states are not orthogonal, {Bh\Bl) ^ 0, which could not occur in a classical system. 
The theoretical prediction of \q/p\ requires the calculation of inclusive nonleptonic rates, which 
is feasible in the heavy-quark limit, but in practice the uncertainties may be sizable. This type 
of CP violation is well established in the charge asymmetry of semileptonic K'^ decay (given by 
6l ~ 2ReeR-) or the asymmetry between K^^q — > e^X and K^^q —>■ e~ X (given by4Reei^) 0. 
The latter asymmetry measuring CP violation in mixing is denoted by A^h in the B system, 
where it is expected to be below the 10~^ level. 

When both B^ and B^ can decay to the same final state, /, the time-dependent CP asymmetry 
can be studied: 

Af{t) = ^ /] - ^[-^°(^) ^ -^j = Sf sin(Amt) - Cf cos(Amt) , (24) 

' T[B^{t) ^ f]+T[B^{t) ^ f] ^ ^ I } \ 1^ \ I 
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where Am = mB^ — tubi^ and t is the proper decay time of the and where we have assumed 
CPT invariance and neglected hfetime differences in the neutral i?-meson system (see Refs. 



and 1 18 1 for a more general formulation of the time dependence). The S and C coefficients can 
be calculated in terms of the B^B^ mixing parameters and the decay amplitudes, 

_ 2ImA/ _ 1-IA/p ^ - l^f 

It is possible that \q/p\ = 1 and |A/| = 1, i.e., there is no CP violation in either mixing or 
decay, but the CP asymmetry in Eq. (]2^ ) is nonzero, because ImAj 7^ {CP violation in the 
interference between decay with and without mixing, or mixing-induced CP violation). 

For decays in which the final state is a CP eigenstate and amplitudes with one weak phase dom- 
inate, Aj^p measures a phase in the Lagrangian theoretically cleanly, independent of hadronic 
matrix elements. (However, / does not have to be a CP eigenstate for this type of CP violation 
to occur, nor for it to have a clean interpretation.) For such decays, the CP symmetry of QCD 
implies Af^p = r]f^pAj^^, where r]f^p is the CP eigenvalue of /, simplifying the evaluation of 
Xf. In such cases Cj^p = and Sf^p = ImA/^p = sin(arg A/^p), where argA/^p is the phase 
difference between the B^ — > / and B^ — > B^ f decay paths. For the generic case, where sev- 
eral amplitudes with different weak phases contribute, detailed knowledge of Af /Af is necessary 
to interpret the experimental measurements. In particular, the strong phases of the amplitudes, 
which usually result in hadronic uncertainties in the interpretation of 5j, and may also give rise 
to 7^ 0, must be determined. Note that if is small, it does not imply that Sf provides 
clean information on short-distance physics. If there are amplitudes with different weak but 
small relative strong phases, then C/ ~ 0, but Sf still depends on the hadronic physics. 



1.5 Physics Beyond the Standard Model 

There are many reasons to believe that there is physics beyond the SM. The evidence for dark 
matter implies not-yet-seen particles, the observed matter-antimatter asymmetry of the Universe 
implies not-yet-observed CP violation, and neutrino masses make the existence of heavy right- 
handed neutrinos an appealing scenario. More aesthetic reasons include the gauge hierarchy 
problem that requires new physics to stabilize the huge Mpianck/^^weak ratio, grand unification of 
the electroweak- and strong-interaction couplings (that appears naturally with Supersymmetry 
(SUSY) at the TeV scale), the strong CP problem that may imply the existence of an axion, 
and the evidence for dark energy that may simply be a cosmological constant in the Einstein 
equations. Some of these may also relate to new flavor physics that could have observable effects 
in low-energy experiments. 

It is interesting to probe the flavor sector, because it is not well understood. In the SM, the 
origin of quark-mass and mixing-angle hierarchies are not understood. And if there is new physics 
at the TeV scale, as conjectured from the gauge hierarchy problem and grand unification, we do 
not understand why it does not show up in flavor physics experiments. A four-quark operator 
(sd)'^ / A^p with 0(1) coefficient would give a contribution exceeding the measured value of 
ek unless tInp ^ 10^ TeV. Similarly, (d6)^/yl^p yields Am^ above its measured value unless 
^NP ^ 10'^ TeV. Extensions of the SM typically have many new sources of CP violation and 
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flavor-changing neutral-current interactions, and in most scenarios, to satisfy the experimental 
constraints, specific mechanisms must be introduced to suppress these. 

Supersymmetric models are a good example. Compared with the SM, their flavor sector con- 
tains 59 new CP-conserving parameters and 41 new CP-violating phases |19,2C|, and many of 
them have to be suppressed to not contradict the experimental data. The smallness of ek and 
AniK have driven SUSY model building to a large extent, and have given rise to several mecha- 
nisms that constrain the relevant soft SUSY-breaking terms. Generic models with superpartner 
masses at the TeV scale also violate the current bounds on electric dipole moments, and fine 
tuning at the 1% level is required to avoid conflict with the data. 

We do not know whether there is new flavor physics at the TeV scale, so it is hard, if not 
impossible, to make firm predictions. If the new flavor scale is at the TeV scale (where the 
hierarchy problem suggests that some new physics should exist), it could have observable effects 
in the (low-energy) flavor sector. If the flavor scale is much above the TeV scale, as is the case for 
example in models with gauge-mediated SUSY breaking, then we expect no signiflcant deviations 
from the CKM picture. In any case, if we see new particles at the Large Hadron Collider (LHC), 
it will be much clearer how to use flavor physics to learn about the new couplings. 



2 Theoretical Framework and Tools 



With the remarkable exception of the UT angles, the experimental observables presently used to 
constrain p and rj depend on hadronic matrix elements. QCD is well established as the theory of 
strong interaction, and it has been tested to high precision in the perturbative regime where the 
coupling constant as is small. However, presently it is difficult to obtain quantitative predictions 
in the low-energy regime, except for a few special cases. Although it is beyond the scope of this 
review to discuss the wealth of approaches to nonperturbative QCD, it is useful to recall a few 
general techniques to evaluate the matrix elements relevant to quark-flavor physics. 

The methods reviewed below all give controllable systematic errors, by which we mean that the 
uncertainties can be incrementally improved in a well-defined way, expanding in small parameters 
order by order. Most of the model-independent theoretical tools utilize that some quark masses 
are smaller while others are greater than vIqcd (throughout this review TIqcd denotes a typical 
hadronic scale, of order 500 MeV). Expanding in the resulting small ratios can be used to simplify 
some of the hadronic physics. However, depending on the process under consideration, the 
relevant hadronic scale may or may not be much smaller than mi,. For example, /^r, mp, and 
m\/ms are all of order tIqcd; but their numerical values span an order of magnitude. In most 
cases experimental guidance is necessary to determine how well the expansions work. 

2.1 Effective Hamiltonians for Weak Decays 

All flavor-changing interactions (except that of the top quark) are due to tree and loop dia- 
grams involving heavy virtual particles: W bosons in the SM, or not-yet-discovered particles 
in its extensions. These particles propagate over much shorter distances than 1/mh, so their 
interactions can be described by local operators. In principle, there is an inflnite number of 
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such operators. The contributions of the higher dimensional ones are however suppressed by 
increasing powers of mf,/mw, so it is sufficient to consider the first few operators. The effective 
weak Hamiltonian can be written as Hw = J2 Ci{fj.) Oi{fi), where Oi are the lowest dimensional 
operators contributing to a certain process and Cj are their Wilson coefficients, with perturba- 
tively calculable scale dependences. It is possible to sum the large logarithms of m'^r/ml by 
first calculating Cj(/i ~ mw) (matching) and then evolving the effective Hamiltonian down to a 
scale /U ~ mf, using the renormalization group (running). This shifts the large logarithms from 
the matrix elements of Oi to Cj. In this last step the operators can in general mix; for example, 
C2{mw) affects all of Ci_9(mfc) below. 

The simplest examples are semileptonic decays, where integrating out the virtual W m. b ^ 
clU, for instance, gives 



w 



Gf 
V2 



{cb)v-A [ly) 



V-A , 



(26) 



where V — A denotes the Dirac structure 7^(1 — 75) between the fermion fields. In this case there 
is only one operator, and its coefficient is scale independent because the (axial) vector current is 
(partially) conserved. Semileptonic decays involving a ii pair and nonleptonic decays are more 
complicated. The AB = 1 effective Hamiltonian has both AS* = and 1 terms. The AS = 1 
part, 

Gf 
V2 



Hw 



E 

p=u,c 



(27) 



j>3 



contains the operators 



Of 


= {uf3ba)v~A{SaUi3)v-A, 


01 


= {ub)y^A {su)y-A , 


01 


= {ci3ba)v-AiSaCf3)v-A, 




= (cb)y.A {sc)y-A , 


03 


= T,q(.sb)v-A {qq)v-A, 


O4 


= Eq(S/3^a)y-A {qaqp)y-A, 


05 


= llq{sh)v-A {qq)v+A, 


Oe 


= Y.q{Sf3ba)v-A {qaqp)v+A , 


07 


= -gf^mftsa'^'^F^, (1+75)6, 


Os 


= -^m,-sa^^-Gl,T-{l + ^^)b, 


09 


= ^{sb)v-Amv, 


do 


= ^{sb)y^A{U)A. 



(28) 



The AS* = part is obtained by replacing s ^ d in Eqs. ( p^ ) and 
current-current operators, O3-6 are penguin operators with a sum over q 



Here 05*2 and Of 2 



are 

u,d,s,c,b flavors 

(a, /3 are color indices). The Ci coefficients in Eq. (27) are known at next-to-leading logarithmic 
order In nonleptonic B decays, four more operators corresponding electroweak penguin 

diagrams contribute to Hw- 



Of 

^9 



Eg 1 eq {sb)y^A {qq)v+A , 

T,q i eg {sb)y-A {qq)v-A , 



^10 



2 



q 2 1 



eq {Sl3ba)v-A {.qaqi3)v+A , 
Si3ba)v-A {qaqi3)v-A ■ 



(29) 



Non-SM physics can either modify the coefficients Ci at the weak scale, or give rise to ad- 
ditional operators (for example, those obtained by replacing (s6)y_A by {sb)y+A)- The goal 
is to find out if the data show evidence for either type of modification. The largest source of 
theoretical uncertainty usually stems from our limited ability to compute the hadronic matrix 
elements of such operators, and the rest of this section introduces methods that allow calculating 
them, or relating them between various observables. 
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2.2 Chiral Symmetry 

The n, d and s-quark masses are small compared with vIqcDj so it is useful to consider the nig 
limit {q = u, d, s) and treat corrections perturbatively. This is known as the chiral limit, because 
the Lagrangian for the light quarks has a SU(3)i x SU(3)ij chiral symmetry, under which the 
left- and right-handed quarks transform differently. This symmetry is spontaneously broken 
to SU(3)v/ by the vacuum expectation value of the quark bilinears, {(fjiQi) = vS'^K The eight 
broken generators are related to Goldstone bosons, the three pions, four kaons, and the r/ (for 
now, we neglect rj-rj'-ir^ mixing). Chiral symmetry relates different hadronic matrix elements to 
one another, and has very diverse applications in flavor physics. 

Because the u and d-quark masses are small, the SU(2) isospin symmetry between the u and 
d is usually a very good approximation. The corrections to the chiral limit are suppressed by 
{nid — m„)/yl;^sB — where vl^gB ~ 1 GeV (of order nip or 47r/7r) is the chiral symmetry breaking 
scale — and are usually not larger than a few percent at the amplitude level. There are also 
explicit violations of chiral symmetry, for example, due to weak or electromagnetic interactions. 
In certain cases these effects can be enhanced as is the case, e.g., in the neutral versus charged B 
(or D) meson lifetimes. The full SU(3) symmetry is broken by mg/yl^^sB) and is known to have 
typically 20 — 30% corrections. (The same is true for its u-spin and d-spin SU(2) subgroups, 
which act on the ds and us pairs, respectively.) 

Some of the most prominent cases of isospin symmetry in the context of the CKM matrix 
include relations between amplitudes involving charged and neutral pions, the determination of 
iKidI (Sec. p. ID , and the extraction of the UT angle a from B ^ tttt decays (Sec. ^). Simi- 
larly, SU(3) symmetry and chiral perturbation theory are key ingredients in determining \ Vus\ 
(Sec. p.2.2| ). It also relates form factors and certain matrix elements involving pions and kaons 
to one another, a relation that has many applications. Recently, the SU(3) relations between 
nonleptonic decays have been extensively studied, because the tttt, Ktt, and KK amplitude re- 
lations give sensitivity to the UT angle 7 and possibly to new physics. SU(3) has also been used 
as a bound on the SM-induced deviations of the time-dependent CP asymmetries from sin 2/3 



in the penguin-dominated modes (see Sec. 7.2) 



2.3 Heavy-Quark Symmetry and Heavy-Quark Effective Theory 

In mesons composed of a heavy quark and a light antiquark (plus gluons and qq pairs), the 
energy scale of strong interactions is small compared with the heavy-quark mass. The heavy 
quark acts as a static point-like color source with fixed four-velocity, which cannot be altered 
by the soft gluons responsible for confinement. Hence the configuration of the light degrees of 
freedom (the so-called "brown muck" ) becomes independent of the spin and flavor (mass) of the 
heavy quark, which, for Nf heavy-quark flavors, results in a SU(2A'^^) heavy-quark spin-flavor 



symmetry [22|. 



Heavy-quark spin-flavor symmetry has many important implications for the spectroscopy and 
strong decays of B and D mesons (for a review, see e.g. |^^). It is especially predictive for 
exclusive B D^*^ii7 semileptonic decays, which are relevant for the determination of \Vcb\ 
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(Sec. 3.2.4). When the weak current suddenly changes the flavor (on a time scale <C Aq 



^QCD^ 



momentum, and possibly the spin of the 6-quark to a c-quark, the brown muck only notices that 
the four-velocity of the static color source has changed, Vb — > fc- Therefore, the form factors that 
describe the wave-function overlap between the initial and final mesons become independent of 
the Dirac structure of weak current, and depend only on a scalar quantity, w = Vf,-Vc- Thus all six 
B D^*^itV form factors are related to a single Isgur-Wise function, ^{vh-Vc), which contains all 
the low-energy nonperturbative hadronic physics relevant for these decays. Moreover, ^(1) = 1 
because at zero recoil — where the c quark is at rest in the b quark's rest frame — the configuration 
of the brown muck does not change at all. 

Deviations from the heavy-quark limit can be included using the heavy-quark effective theory 
(HQET) [^4|, which provides a systematic expansion in powers of as(mQ) and ^qcd/'"^q {Q = 
b,c). The former type of corrections are calculable perturbatively, whereas the latter ones can 
be parameterized by a minimal set of hadronic matrix elements that can be extracted from 
data and/or estimated using nonperturbative techniques. Heavy-quark spin symmetry also 
implies relations, in combination with chiral symmetry, between, for example, B pil7 and 
B K*i+e- form factors |5|. 

2.4 Factorization and Soft-Collinear Effective Theory 

Researchers have long known that in the decay B M1M2, if the meson Mi that inherits 
the spectator quark from the B is heavy and M2 is light then "color transparency" can justify 



factorization |26, 27 1. Traditionally, naive factorization refers to the hypothesis that matrix 
elements of the four-quark operators can be estimated by grouping the quark fields into a pair 
that can mediate B Mi transition and into another pair that describes vacuum M2 
transition. 

Recently the development of the soft collinear effective theory (SCET) | p8| , |2^ put these 
notions on a firmer footing. SCET is designed to describe the interactions of energetic and low 
invariant-mass partons in the Q ^ ^qcd limit. It introduces distinct fields for the relevant 
degrees of freedom, and a power-counting parameter A. There are two distinct theories, SCETj 



in which A = ^Aqcd/Q and SCETn in which A = Aqqd/Q. They are appropriate for final 
states with invariant mass QX; i.e., SCETj for jets and inclusive B — > Xg'j, Xuii', Xs£~^£~ decays 
(m^ ~ AqcdQ), and SCETn for exclusive hadronic final states {m? ~ ^qcd). It is convenient 
to use light-cone coordinates, decomposing momenta as 

= {n-p)—+p^^ + {n-p)— =p-—+p± ('^^) 

where = = and n-n = 2. For a light quark moving near the n direction, p_ ^ pj_ ^» p+, 
and this hierarchy can be used to successively integrate out the less and less off-shell momenta. 
Factorization has been proven in {B D7r^)-type decays to all orders in a<j and at leading 
order in Aqcd/Q by decoupling the ultrasoft [{p^,p^,p_i_) ~ Q{X^ , X'^ , X'^)] gluons from the 
collinear [{p-,p+,p±) ~ Q(l, A^, A)] Lagrangian at leading order in A. 

The study of heavy to light transitions, which are particularly important for CP violation and 
CKM measurements, is complicated by the fact that one has to work to subleading order in A. 
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This is because collinear and ultrasoft quarks cannot interact at leading order, only via the mixed 
ultrasoft-collinear Lagrangian, C^^J , which is suppressed by one power of A and allowed to couple 



an ultrasoft and a collinear quark to a collinear gluon |31|. This is relevant for all processes 
in which the spectator quark in the B ends up in an energetic light meson. So far, maybe the 
most surprising model-independent result at subleading order is proof of factorization in the 
previously intractable B D^ir^ type "color-suppressed" decays [^2|. In this case factorization 
means the systematic separation of the physics associated with different momentum scales, and 
the factorization theorem can even accommodate a nonperturbative strong phase. 

Other important applications include the exclusive semileptonic form factors, and factorization 
in charmless two-body decays, both of which are subject to intense discussions. For semileptonic 
B Tidy and ptu decays, there are two contributions of the same order in ^qcd/Q- It is not yet 



clear whether or not the factorizable parts that violate the form factor relations |33| are smaller 



than the nonfactorizable contributions. For charmless decays [|34| , p5|| , an additional complication 



is that the power counting appropriate for the charm penguins has not been settled [|36| , |37| . 
2.5 The Operator Product Expansion for Inclusive Decays 

In the large mf, limit, there is a simple argument based on a separation of scales, that inclusive 
rates may be modeled by the decay of a free 6-quark. The weak decay takes place on a time 
scale much shorter than the time it takes the quarks in the final state to hadronize. Once the 
6-quark has decayed, the probability that the decay products hadronize is equal to unity, and 
the (unknown) probabilities of hadronization to specific final states can be ignored. 

This argument can be formalized using an operator product expansion (OPE) |3^. To cal- 
culate the semileptonic rates, the leptonic and the hadronic parts of the {Xiv\H \ B) matrix 
element can be separated. Nonperturbative strong-interaction effects enter the rate via the 
hadronic tensor, 

;^.. = ^p.)=,.(,,_,_„,M>(™, (31) 

where J'^ = q'j^Pib for semileptonic B decay. The optical theorem implies that W^'^ can 
be related to the discontinuity across the cut of the forward-scattering matrix element of the 
time-ordered product, 

J 2mB ^ ^ 

i.e., W'^'^ = —^ImT^'^. When the energy release to the hadronic final state is large, the time- 
ordered product can be expanded in local operators, whose i3-meson matrix elements can be 
parameterized using HQET. T^^'^ is an analytic function with singularities that correspond to 
on-shell hadronic intermediate states. For semileptonic decays one can deform the integration 
contour away from the physical region in the complex q ■ v plane (see Fig. P), and asymptotic 
freedom ensures that T^'^ can be reliably calculated far from its singularities. 

In the rrify ^ ^qcd limit, the OPE reproduces the 6-quark decay result, and shows that the 
leading power-suppressed corrections occur at order ylq^j-j/m^ and can be parameterized by two 
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Re q • V 



Figure 2: (Left) Forward-scattering amplitude whose imaginary part gives the decay rate. (Right) 
Analytic structure of the correlation function in the complex q-v plane and the integration contour 
that gives the total rate. 



matrix elements, usually denoted by Ai and A2. The result is of the form 

0(Ai,A2) 
1 + — + ^ +■■■ 

nib nT'b 



^-p _ / 6-quark \ 
decay J 



(33) 



Most decay rates of interest have been computed including the ^qcd/™'6 corrections (which 
are parameterized by six additional hadronic matrix elements), whereas the 6-quark decay rates 
have been computed including perturbative corrections at order and ct^/^o (/3o = 11 — 2nj/3 
is the first coefficient of the QCD /3-function, and the terms proportional to it often dominate 
at order a^). 

The OPE can also be used to calculate inclusive nonleptonic decays, such as heavy hadron 
lifetimes, the Bg width difference, Ar^, and the semileptonic CP asymmetry, ^SL- Whereas the 
structure of the results is the same as for semileptonic decays, there is no "external" current 

to the strong interaction in these cases, so there is no variable that can be used to move the 
integration contour away from the cuts. Therefore, the OPE is performed directly in the physical 
region, and we expect the accuracy of these calculations to be worse than those for semileptonic 
decays. Experimental data are crucial to validate each of the applications; in particular, the 
unexpected smallness of the Ab to B lifetime ratio may indicate that the nonleptonic calculations 
are more difficult to control. 



2.6 Lattice QCD 

In lattice QCD (LQCD), the functional integral for correlation functions is computed using nu- 
merical Monte Carlo integration. To keep the degrees of freedom and the required computational 

power at an acceptable level, a finite lattice spacing, a, is introduced to discretize space-time, 
and the calculation is performed within a finite volume. The finite lattice spacing (finite volume) 
is an ultraviolet (infrared) cutoff, which can be dealt with using field theory techniques to extract 
cutoff-independent quantities from the results of the computations. In practical calculations the 
ultraviolet cutoff is typically approximately ^ (2-3) GeV. Consequently, 6-quarks can be 
simulated only using actions based on HQET and nonrelativistic QCD, where rub is removed 
from the dynamics using the effective theory. The matching between the lattice and the contin- 
uum is in principle straightforward, but in practice it often introduces a significant uncertainty, 
because as{a~^) is not small, and few matching calculations have been done beyond one loop. 
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The most difficult issues are related to how the fermionic part of the action is computed. Until 
a few years ago, to make the computations manageable, most calculations used the quenched 
approximation, which amounts to neglecting virtual quark loops. It is not a systematic ap- 
proximation to QCD and the error associated with it is difficult to estimate. Recent advances 
in algorithms and computers have allowed an increasing number of quantities to be computed 
with dynamical u, d, and s-quarks. A remaining source of systematic uncertainty is often the 
chiral extrapolation. The simulations with light u and d-quarks are computationally expensive, 
so the use of m^^d/ms ~ 1/2 is not uncommon. To date only the MILC collaboration has used 
niu^d/iT^s ~ 1/8, which still has to be extrapolated to the physical values, mu^d/iT^s ~ 1/25. The 
results based on MILC's gauge configurations use improved "staggered" quarks. There are sev- 
eral ways to discretize the quark action, and the cost of simulations with light quarks is sensitive 
to this choice. Unquenched calculations with staggered quarks are much faster than with the 
Wilson, domain wall, or overlap formulations. In the staggered formulation each field produces 
four quark "tastes", so there are 16 light mesons. The unphysical ones are removed by taking 
the fourth root of the fermion determinant. The validity of this fourth-root trick has so far been 
proven only in perturbation theory. 

Recently many new results have emerged in which all these systematic effects are studied 
in detail. The gold-plated quantities are those that contain at most one hadron in both the 
initial and final state, not moving with large velocities. These include decay constants, bag 
parameters, and scmileptonic form factors, which are important for flavor physics. For most of 
these quantities, the currently available calculations using staggered fermions obtain significantly 
smaller errors than those with Wilson and other formulations. The fastest LQCD computers 
now run at a computation rate of a few tens of TeraFLOPs, and the PetaFLOP barrier is 
expected to be passed by the end of this decade. This may allow LQCD to make predictions 
approaching the percent level, especially if all ingredients of the calculation can be performed 
numerically (including matching from the lattice to the continuum). With such precision, several 
measurements in exclusive decays, which are sensitive to new physics but are presently dominated 
by theoretical uncertainties in the corresponding matrix elements, could become precise tests of 
the SM (such as B pj, B — > ir£i>, £k, etc.). 

3 Determining Magnitudes of CKM Elements 

Determining of the magnitudes of CKM elements uses a number of sophisticated theoretical and 
experimental techniques, the complete discussion of which is beyond the scope of this review. 
We concentrate on topics where recent developments have occurred. 

3.1 \Vud\ from /3-Decays 

The CKM matrix element \Vud\ has been extracted by three different methods: superallowed 
nuclear /3-decays with pure Fermi transitions (O"*" O"*"), neutron /3-decay (n peUf.), and 
pion /3-decay (vr"'' t:^ e'^ Veij)) ■ All three methods involve the hadronic form factor of the 
vector current (/|li7^(i|i), where |i) and |/) are the initial and final states of the transitions. 
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The neutron /3-decay also involves an axial- vector form factor, which requires external input for 
the ratio of the axial-vector to the vector coupling constants. The form factors describe the 
long-distance hadronic effects that confine the quarks within hadrons. The normalization of 
the vector current is fixed at the kinematic endpoint of zero-momentum transfer between initial 
and final states, and at the exact isospin limit {rfXu 

= md). Isospin-breaking corrections are 
suppressed by the quark mass difference divided by the hadronic scale. 

Superallowed nuclear /3-decays give the best experimental precision. The product of the 
integral over the emitted electron energy spectrum and the neutron lifetime is proportional 
to ll^dl"^. Radiative and charge symmetry breaking corrections, which depend in part on the 
nuclear structure of the nucleus under consideration must be applied. They lie between 3.1% and 
3.6% for the nine superallowed transitions that give the best sensitivity |3^. Recent theoretical 
improvements in the calculation of loop contributions, notably for the problematic box 
diagram, and a detailed error analysis achieved a reduction in uncertainty by a factor of two to 
1.9 X 10"^ m. With this, the world average is |Kd|o+^o+ = 0.97377 ± 0.00027 ||. 

From a theoretical point of view, the neutron /3-decay, or ultimately the pion /3-decay, represent 
the cleanest determination of \ Vud\-, because they are free from nuclear structure effects. Because 
no other decay channels exist, the decay rate can be obtained from lifetime measurements. The 
PERKEO-II experiment achieved significant improvement in this \Vud\ determination owing to 
the precise measurement of qa/qv = —1-2739 ib 0.0019 Q| (confirmed in Ref. [^]), a value 
that is, however, significantly larger than earlier measurements. Several new experiments are 
planned or are in construction to remeasure this quantity. Unfortunately, there is currently 
a problem with the second input to determine \Vud\ from neutron /3-decays, namely a large 
discrepancy between a new neutron lifetime measurement [E3] and the PDG average |45|. Such 



a discrepancy, if confirmed, would significantly impact the extracted value for \Vud\^ and hence 
CKM unitarity. Ignoring this measurement, the present world average from neutron /3-decays is 



\Vud\nfi = 0.9730±0.0004±0.0012±0.0002 gl[, where the errors arise from the neutron lifetime, 
from gA/gVi and from radiative corrections, respectively. Note that the change induced by the 
new lifetime measurement, if confirmed, would be as large as three times the total error quoted. 

The pion /?-decay Tr"*" — > -K^e^Uf, is an attractive candidate to extract \Vud\ from the branching 
ratio of the decay and the pion lifetime. It is mediated by a pure vector transition. Unfortunately, 
because of the small branching ratio of 10~*, the statistical precision is not yet competitive with 
the other methods: \Vud\ = 0.9748 ± 0.0025 p|,p|. 



3.2 CKM Elements from (Semi)Leptonic Decays 
3.2.1 Leptonic Decays 

The decay rate of the lightest pseudoscalar meson, M — which is composed of the quarks qdQu — to 
a lepton antineutrino pair is given by 

r (M- ^ rv,) = ^ \V,^JflmMm} [l - , (34) 

where the decay constant /m, defined by ifhip^ = (0|gu7^75grf|M(p)), contains all the nonper- 

turbative strong-interaction physics. If /m is known, then leptonic vr, D, Dg, B, and Be 
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decay rates determine iKdl, \Vcs\, |Kb|, and \Vcb\, respectively. 

Owing to the mj suppression, however, these rates are small for £ = fi or e, and the r final 
state is harder to reconstruct experimentally. The current errors of the LQCD determinations of 
the decay constants are around the 10% level, so the observed leptonic K, D, Dg and B decays 
do not give competitive direct determinations of the corresponding CKM elements. However, 
ratios of decay constants are easier to determine in LQCD, and a recent calculation of JkI^-k = 
1.198 lb 0.0031qqq5 [47| allows |Kis/Ktd| to be extracted from the ratio of leptonic K and vr 



decays, yielding |Ks| = 0.2245l°:[J[]^2 |8|. The ratios fosl fo and are also known much 

more precisely than the individual decay constants. However, the experimental errors are still 
sizable for D — > £v, and because the Bg meson is neutral, it cannot decay to a tP pair.^ 

Metrologically useful experimental information exists for the decay i?^ — > t~^i^t, where av- 
eraging the experimental likelihoods from Belle (reporting a 4.2(7 evidence for this decay) and 
BABAR gives (10.4 Ig'y) ^ 10~^ for the branching fraction, which is consistent with the pre- 
diction, (9.6 lb 1.5) X 10~^, from the global CKM fit (not including the direct measurement) 
It is interesting to study the constraint obtained from this measurement in conjunction with 
the measurement of the B^B^ oscillation frequency, Am^, which also depends on fs, so that 
fs cancels in the combination. The allowed regions obtained in the (p, rj) plane are shown in 
Fig. ^. The remaining theoretical uncertainty stems from the bag parameter, Bb, that enters 
the SM prediction of Am^ [see Eq. (li! 



3.2.2 

Traditionally, the magnitude of Vus has been extracted from semileptonic K decays, which can 
be analyzed using chiral perturbation theory. The K — > Triv amplitude can be expanded in 
powers of the K and vr momenta and the quark masses divided by a scale of order ^^SB- The 
decay rates for K = K^, and i = e, ^ can be written as 



T{K ^ nev) = ^g^CI^' |/f ^10)|' (1 + . . .) , (35) 

where C = 1 (1/2) for (K^) decays, I^^ is a phase-space integral depending on the shape of 
the form factor, f^{q'^), and the dots stand for isospin-breaking and electroweak corrections. It 
is convenient to factor out f^"'" (0) = /+(0), because the Ademolo-Gatto theorem |5^ implies 
that there is no 0{p^) correction to /+(0), i.e., /+(0) = 1 -b /p4 -b fp6. Most of the ©(p^^) 
corrections are known in terms of physical parameters, and beyond this the estimates are model 
dependent. The original value, /+(0) = 0.961 ib 0.008 is still broadly accepted, although 
some other calculations differ by as much as 2%, with uncertainties around 1% |4^. It will 
be interesting to see whether unquenched LQCD results confirm this estimate. Significant 
progress on the experimental side has taken place recently. After a high-statistics measurement 
of the B{K^ — > TT^e^u) a suite of measurements of neutral-kaon branching ratios 

^The hadronic physics relevant to the rare decay Bs is also governed by /s^ , and once the decay is 

measured at the LHC, it will give a precise determination of \Vub\ from the double ratio [B{B £V)/B{Bs 
£+£~)] X [B{Ds eu)/B{D ^ eu)] |49|. This double ratio can also give a precise SM prediction for Bs l+T, 
and B ^ tU can even be replaced by _B — > t^t~ . 
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Figure 3: Constraints of 95% C.L. in the {'p,r]) plane from the t~^Vt- branching fraction, 

Amrf, and the combined use of these j5i| / (darker shades indicate 68% C.L. regions) ]5^. The 
results are compared with the global CKM fit (cf. Sec. |^. 



of form factors [56|, and of lifetime [57| followed. From semileptonic kaon decays one obtains 



IK. 



0.2257 ±0.0021. 



The T spectral function into jAS"! = 1 final states can be used to determine the s-quark mass 
and the CKM element \Vus\ |5^- Although the determination of rns{m.j-) is strongly model 
dependent, the determination of \Vus\ appears on safe grounds, because it relies on the non- 
weighted spectral function with maximum inclusiveness and a better control of the perturbative 
expansion. Present data give || |T4s| = 0.2204 ±0.0028 ±0.0003 ±0.0001, where the first error 
is experimental, the second theoretical and the third due to (taken from LQCD). 



3.2.3 I Fes I and \V^\ 

For the dominant semileptonic D decays, D K^*^ii', neither heavy-quark nor chiral symmetry 
is particularly useful. Hence, until recently, the only precise direct constraint involving \ Vcs\ was 
from the leptonic branching ratio of the W, which measured the sum of CKM elements in the first 
two rows (corrected for effects from perturbative QCD), J2i,je{u,c,d,s,b} — l-999±0.026 
(compared with the expected value of two), providing a precise test of unitarity. 

Fortunately, the D Kiu and Triv decays are well suited for LQCD calculations, because 
the charm quark is neither too heavy nor too light, and the maximal momenta of the final state 
X or TT is not too large. Together with the CLEO-c program, it promises significant advances 
both in the measurements of \ Vcs\ and \ Vcd\ as well as in our understanding of the robustness of 
LQCD calculations from which B physics will also benefit. 

The direct determination of \Vcs\ and \Vcd\ is possible from semileptonic D or leptonic Dg 
decays, relying on the form factors (which depend on the invariant mass of the lepton pair, g^). 
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Recently, the first three-flavor unquenched LQCD calculations for D KlV and D ttIv were 



published |61|. Using these calculations, and the isospin-averaged charm semileptonic widths 



measured by CLEO-c, one obtains [£2| 

iVcdl = 0.213 ±0.008 ±0.021, jVcsl = 0.957 ± 0.017 ± 0.093 , (36) 

where the first errors are experimental, and the second theoretical due to the form factors.^ 
Note that both \ Vcs\ and \ Vcd\ are determined much more precisely if CKM unitarity is assumed, 
with uncertainties around 0.0002 and 0.001, respectively (see Table |3| on p. |43| ). Therefore, these 
measurements do not affect the global CKM fit at present. 

3.2.4 I Feb I 

The CKM element \Vcb\ can be extracted from exclusive and inclusive semileptonic B decays. 
The determination from exclusive B D^*^itV uses an extrapolation of the measured rates to 
zero recoil, w = 1 (i.e., maximal momentum transfer to the leptons, Qmax)- The rates can be 
written schematically as 

dr(S^L>Wa7) ^ , ,^^,2 [(u;2-l)i/2^2(^) ioi B^D*, 

— ^ = (known terms) x V^ft >^ { , ^ Lio o) [ (37) 



Heavy-quark symmetry (see Sec. p.3[ ) determines the rate w = 1, that is, in the mq » ^qcd 
limit {Q = 6, c), J^{w) = J^^^iw) = £,{w), and in particular .F(^,)(l) = 1, allowing for a model- 
independent determination of \Vcb\- The corrections to the heavy-quark limit can be organized 
in expansions in as{mQ) and ylQCD/'^T-Q, which take the form 

^ , , , , (LQCD or models) 

= 1 ± CA as ± ± ^^^^ 2 + • • • ' 

^, . , , (LQCD or models) , , 

^(1) = l + cy(Q,)±^-:^ ^±.... (38) 

mq 

Heavy-quark symmetry implies that the leading-order results are unity, whereas the absence of 
the ^qcd/w-q correction to T^iX) is due to Luke's theorem [Q. The perturbative corrections 



CA = —0.04 and cy = 0.02 are known to order oq, |64|, and higher-order terms are below the 1% 
level. The terms indicated by "LQCD or models" are known only from quenched LQCD ||6^ or 
phenomenological models. The experimental results for the rates, extrapolated to if = 1 are 



|ycb|.F*(l) = (37.6 ±0.9) X 10"^ |ycfe|^(l) = (42.2 ±3.7) X 10"^ (39) 

Using T^iX) = 0.91 ± 0.04, yields \Vd\ = (41.3 ± l.Ocxp ± 1.8th) x 10~3_ xhe B DtU data 
are consistent with this, but to make a real test, smaller experimental errors and unquenched 
LQCD calculations are needed. 



The inclusive determination of \Vcb\ is based on the OPE (see Sec. 2^). The state of the 



art is that the semileptonic rate, as well as moments of the lepton energy and the hadronic 

^The most precise measurement of IVcdl still comes from charm production in neutrino scattering. The PDG 
quotes |Vcd| = 0.230 ± 0.011 an uncertainty that will be hard to reduce, whereas that from Eq. ( ^6|) will 
improve with better LQCD calculations and data. 
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invariant-mass spectra with varying cuts have been computed to orders ^qcd/ n^b and a^Po. 
Their dependence on m^^c and the parameters that occur at subleading orders in ylQCD/"^6 are 
different. The expressions for a large number of observables are of the form 



T(B ^ X^eu) = (^] (0.534) [l - 0.22 ( 

^ ^ 1927r3 \ 2 J ^ ^[ V500 MeV / 



0.011 ( ^J^]' - 0.052 ( .^J^] - 0.071 ( 



SOOMeVy U500 MeV)2y ^500 MeV) 



0.096. - 0.0304lm + 0-015ef ^^^^ ) + • • • 



(40) 



where mr is the T{1S) mass. The precise determination of \Vcb\ requires the use of an appropriate 
short-distance quark mass, the so-called "threshold" mass |6^,^, and Ais = mr/2 — m]f is 
related to one of these, [37,^|. The ylq^j-j/mf corrections are parameterized by Ai,2. 



The ylq^j-j/m^ terms are known and are parameterized by six more nonperturbative matrix 
elements. For the perturbative corrections e = 1 shows the order in the expansion, and the 
BLM subscript shows that only the terms with the highest power of /?o are known |7^. Such 
formulas are fit to approximately 90 (correlated) observables obtained from the lepton energy 
and hadronic invariant-mass spectra. The fits determine \ Vcb\ and the hadronic parameters, and 
their consistency provides a powerful test of the theory. The fits have been performed in several 
schemes and give ]7l|,^, 

iKfel = (41.7 ±0.7) X 10-3, (41) 

where the central values were averaged and the errors quoted in each paper were retained. The 
same fits simultaneously determine the quark masses, m^^ = (4.68 it 0.03) GeV and rnc{fnc) = 
(1.22 ± 0.06) GeV [0,^], which correspond to mb{rnb) = (4.18 ± 0.04) GeV and ml^ = (1.41 ± 
0.05) GeV. The rub — rric mass difference is even better constrained [fn] , |72| . 



3.2.5 iKbl 

Of the measurements of \Vub\ from exclusive decays, B — > 7r£z7 is the most advanced, as both 
experiment and LQCD calculations are under the best control. The determination relies on 
measuring the rate and calculating the form factor fj^{q^), 

d? ^4^ ^' • 



Unquenched calculations of /+ exist only for > 16 GeV^ (small \ptt\) |74,^, where the 
available statistics is reduced because the phase space is proportional to \'Pt^\^. Averaging the 
LQCD calculations, and using data only in the > 16 GeV^ region, gives \ Vub\ = (4. 13 it 0.62) x 
10-3 pll76|. 



Some of the current \Vub\ determinations use model-dependent parameterizations of f+{(f') 
to extend the LQCD results to a larger part of the phase space, or to combine them with 
QCD sum-rule calculations at small q^ (which tend to give smaller values for \Vub\ |7^)- These 
model-dependent ingredients can be avoided using constraints on the shape of f+{q^) that follow 
from dispersion relations and the knowledge of f+{q^) at a few values of q^ ||78|. The recent 
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LQCD results revitalized this area |7^,|8^. Using the LQCD calculations of /+ at large q'^, the 
experimental measurements, and the dispersion relation to constrain the shape of gives 
\Vub\ = (3.92 ± 0.52) X 10-''^ using the full q"^ range. 

The determination of | Vub \ from inclusive semileptonic B decay is more complicated than that 
of \Vcb\, because of the large B X^tU background. Much of the recent experimental progress 
is related to the fact that the i?-factories produce pure BB pairs, where the full reconstruction 
of one B allows the measurement of the energy and momenta of both the leptonic and the 
hadronic systems in the X^fF decay of the other B (the so-called "S-beam" technique). This 
provides several handles to suppress the B XciV background. The total B —>■ XutP rate is 



known theoretically at the 5% level [37, 38|, but the cuts used in most experimental analyses to 
remove the B X^tU background complicate the theory. The local OPE used to extract \Vcb\ 
is applicable for B X^iV if the typically allowed mx and Ex of the final state satisfy 

m\ > ^xvIqcd » ^QCD • (43) 

However, most of the current \Vub\ determinations restrict the phase space to (a subset of) the 
mx < mn region. As a result, there are three qualitatively different regions. 

If mx ~ ^QCD, the final state is dominated by resonances, and it is not possible to compute 
inclusive quantities reliably. If m\ ~ -Ex^qcd ^ ^qcD' OPE becomes an expansion in 
terms of 6-quark light-cone distribution functions (sometimes termed shape functions). This 
is the case for many measurements to date, such as the rate for > (m^ — m'jj)/{2mB), 
mx < m£i, or Px{= Ex — \px\) < 'm'j^/mB- At leading order one such nonperturbative function 
occurs [^,^], which can be determined from the B Xg^ photon spectrum. At order vlQco/^T-b 



there are several new functions |83], which cannot be extracted from the data, and are modelled. 
The hadronic physics parameterized by functions is a significant complication compared with 
the determination of |V^fe|, where it is encoded in a few hadronic matrix elements. Moreover, at 
order as, moments of the shape function are no longer given simply by local hadronic matrix 



elements i84|, p5||. One can either fit the shape function from B Xg"/, or directly relate 
the B Xu^ partial rates to weighted integrals of the B Xg^ spectrum [^. Finally, if 
m\ » E'x^QCD ^ ^QCD' OPE converges, and the first few terms give reliable results. 
This is the case for the rate for q^ > (m^ — mo)"^ (8^, but the expansion parameter is then 
AqQD/mc- The dependence on the shape function can be kept under control by combining q^ and 
mx cuts. The shape-function dependence can also be reduced by extending the measurements 
into the B XctU region. Recent analyses could measure the B — > X^diy rates for |pe| > 



1.9 GeV and for mx < 2.5 GeV [89|. Averaging these inclusive measurements, the Heavy 



Flavor Averaging Group (HFAG) obtains \Vub\ = (4.38±0.19±0.27) x 10"^ H using [||, which 



indicates a slight tension with the CKM fit for \Vub\ dominated by the sin 2/3 measurement (cf. 
Table I on p. p|). 



3.3 CKM Elements from Rare Loop-Mediated B Decays 

To measure \Vtd\ and \Vts\, one has to rely on loop- mediated rare B decays, or rare K decays 



discussed in Sec. 4.2, or BB oscillations discussed in Sec. 3.4. All of them provide important 
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tests of the SM, as different transitions can receive different contributions from new physics. 

The inclusive branching fraction B{B Xsj) = (3.39 ig-J?) x 10"^ |6) is sensitive to \ VtbVtl\ 
via penguin diagrams with a top quark. A large contribution to the rate stems from the b — > ccs 
four-quark operator, O2 in Eq. ([2^), mixing into the electromagnetic penguin operator, O7. Using 
unitarity, VcbV*, = -VtbVtl-VubV*,, the measured rate implies \VtbVtl\ = (39.0 ±3.1) x 10^^ [||]. 

The theoretical uncertainties can be reduced by taking ratios of processes that are equal in the 
flavor SU(3) limit to determine | V^cj/Vtsl. The exclusive rare decays suffer from larger theoretical 
uncertainties owing to unknown hadronic form factors, but exclusive branching fractions are 
easier to measure. Recently, Belle observed the first b ^ dj signals in i? — > {p/u;)^ decays. We 
define the ratios of branching fractions 



B{B K*j) 



Vt 



td 



Vt 



ts 



' niB - nif, \^ I i(Cv'o^) 2 , for V^-f 



^{ 7 (44) 



rriB - rriK* J \ (^^±7) ^ , for V^-f 



where ^yo^ and Cv^-y denote SU(3)-breaking corrections. The ratio of the neutral rates gives 
the theoretically cleaner determination of |Vtd/Vts|, because weak annihilation enters in the p^'y 
mode, but only W exchange contributes to (Although both annihilation and exchange 

are suppressed by AqcD/nib, they are hard to estimate, and exchange is color-suppressed com- 
pared with annihilation.) Averaging the results from BABAR and Belle |9^,|6^ gives B{B^ — > 
p^-i)/B{B^ i^*°7) = 0.0095 ± 0.0045 |6|. Using the estimate ^yo^ = 1.2 ± 0.1 for this 



average implies 114^/14^1 = 0.16 ± 0.05. Although the theoretical uncertainty for this deter- 
mination of |Vtrf/Vts| will not compete with that for Am^^/Ams, it provides an important test 
of the SM, as new physics could contribute differently to these decays and to BB mixing (cf. 
right plot of Fig. 1). 

3.4 Neutral B-Meson Oscillation 

The B^B^ oscillation frequency is governed by the mass difference Am^ between the two B^ 
mass eigenstates, Bh and Bl. It is defined as a positive number and has been measured by many 
experiments leading to an overall 1% precision, Amrf = (0.507 ± 0.005) ps'^ (dominated by the 
measurements at the i3-factories) [^6|. In analogy to |ex| (see below), B^B^ oscillation in the 
SM is described by FCNC box diagrams. However, in contrast to \£k\^ where the large quark- 
mass-dependent hierarchy in the FCNC amplitude is offset by the tiny CKM matrix element 
|Vf(iV^*|, the IS.B = 2 box diagrams are dominated by intermediate top quarks. This implies that 
Amrf is determined by short-distance physics up to the matrix element of the {d'y'^{l — 75)6)^ 
operator, parameterized by f^^Bd, 

^"^d = ^VBmBafBaBdm^,S{xt) \VtdV*b\'^ . (45) 

Here rjB = 0.551 it 0.007 is a perturbative QCD correction to the Inami-Lim function S{xt) 
(with xt = Wq /m^y) |9|], and Ib^VB^ is taken from LQCD. Much progress has been achieved 



in this domain with unquenched calculations becoming available (see Sec. 2^). A fixed value of 
the CKM factor iVi^V^^p occurring in Eq. (45) describes approximately a circle around (1,0) in 



the {p,r]) plane, to which an elliptic distortion appears only at order A^. 
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In the SM, the mass difference Am^ between the heavy and the hght mass eigenstates 
has only C(A^) dependence on the Wolfenstein parameters p, rj. A measurement of Arrig is 
nevertheless useful for CKM metrology within the SM, because it leads to an improvement in 

tb\ 



the constraint from the Am^ measurement on |Vf^Ktp. Its SM prediction is given by 



^"^^ = ^VBmBA'f'LBdmiySixt) \VtsVtl\' , (46) 



where the parameter ^ = fBs^B^/fB^VBd quantifies SU(3)-breaking corrections to the matrix 
elements, which can be calculated more accurately in LQCD than the matrix elements them- 
selves. Hence a measurement of Am^ reduces the uncertainty of LQCD calculations 
using Wilson fermions have to work with light-quark masses of order 100 MeV, so calculations for 
mesons need to be extrapolated to the chiral (massless) limit. This is not necessary for the 
B^ due to the relatively heavy strange quark. The use of staggered fermions allows LQCD calcu- 
lations to be performed with significantly smaller light-quark masses, which implicates smaller 
systematic uncertainties. 

Because of the large ratio iVts/Vj^P and SU(3)-breaking corrections, B^Bg oscillation occurs 
approximately 35 times faster than B^B^ oscillation. Also, measurements of the time-integrated 
mixing probability Xs indicate a value close to its maximum, 1/2, requiring very fast oscillation. 
Excellent proper-time resolution and large data samples are needed to resolve B^B^ oscillation 
in the detector. Because only pairs of the two lightest B mesons can be produced in T(45) 
decays, B^B^ oscillation can be probed only at high-energy e'^e~ or hadron colliders. 

Lower limits on Am^ have been obtained by the LEP, SLD and Tevatron experiments. A 
convenient approach to averaging the various results is the amplitude method (pSl; see also the 



studies in Refs. ||l3,94]). It consists of introducing an ad hoc amplitude coefficient, A, placed 
in front of the cosine modulation term that describes the time-dependent mixing asymmetry 

( A^unmixod (*)- Admixed (i))/(A^unmixod(i) + Admixed (0) = A-cos{Amst). The advantage of this indi- 
rect probe for oscillation is that the dependence on A is linear so that A is Gaussian distributed. 
The measured amplitude must reach ^ = 1 at the true values of Am<j (if within reach of the 
experimental sensitivity), and zero elsewhere. The preliminary combined 95% C.L. lower limit 
was Anig > 16.6 ps^^ | p6| . 

Shortly before finalizing this review, the D0 and CDF Collaborations reported new prelimi- 
nary results on B^B^ oscillation using data samples corresponding to an integrated luminosity 
of Ifb^^. Using the charge of the muon in a partially reconstructed Bg D~fi'^v^ decay to 
tag the flavor of the B^ meson, D0 sets the 90% C.L. interval 17ps"^ < Am^ < 21 ps"^ [||]. 
Because the upper bound is close to the sensitivity limit of the measurement, larger Am^ values 
are excluded only at the 2.1a level. Using hadronic and semileptonic B^ decays, CDF finds 

Arris 



(l7.3ll0;f8±0.07) ps-i, (47) 



where the first error is statistical and the second is systematic. The probability that the observed 
signal is a background fluctuation is quoted as 0.2% (5%) for CDF (D0). The value of Arris is in 
agreement with the prediction from the global CKM fit, (21.7 142) ps"^, obtained without using 
the measurement in the fit (cf. Table on p. H^). The preliminary world-average amplitude 
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Figure 4: Left: preliminary world- average amplitude spectrum as a function of Ams k6(\], dom- 
inated by the recent CDF measurement j^d^. Right: \Vt(i/Vts\ determined by CDF from the 
ratio Am^/ Arug (full circle) |^/ , from the ratio B{B^ p^'j)/B{B^ — > K*'^'j) (triangle) (see 
Sec. 3.S), and from the global CKM fit (shaded region; not including Arrig in the fit (see Sec. 



scan is shown in the left plot of Fig. ^ The combined signal significance amounts to 3.8a 
at 17.5 ps~^ ||6^. CDF infers from its measurement of Am^ that the ratio of CKM elements 
\Vtd/Vts\ = 0.208 toZl (dominated by the theoretical uncertainty on ^ = 1.21 t[j;{]3^, used by 
CDF). A comparison with the CKM fit is given in the right plot of Fig. ^ This measurement 
provides the first strong indication that B^B^ mixing is probably SM-like. 



4 CKM Matrix Constraints from Kaon Physics 

The neutral-kaon system constrains the UT through K^K^ mixing, through indirect and direct 
CP violation, and through the rare decays — > -k'^vV and the yet unobserved -k^vV. In 



this case, contrary to Eq. (|23D, the convention is to label the states by their lifetimes, \Kg^) = 
Pk\K^) i qK\K^) (experimentally, m^o > w-^o). The SM prediction for Am/< suffers from 
badly controlled long-distance contributions to the mixing amplitudes. Nonperturbative physics 
with large hadronic uncertainties also impede a reliable SM prediction of direct CP violation. 
The above-mentioned rare decays are cleaner and will give precise constraints as soon as they 



are measured with sufficient accuracy. They are discussed in Sec 4.3 



Indirect" CP violation refers to situations in which the CP- violating phase can be chosen to appear only in 
the mixing amplitudes, whereas "direct" CP violation occurs when a CP violating phase must appear in decay 
amplitudes. The observables Imex and Ime^ measure CP violation in the interference between decays with and 
without mixing, whereas ReeK and Ree^, measure CP violation in mixing and in decay, respectively. 
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4.1 Indirect CP Violation 

The most precise measurement of the CP-violation parameter ej^ uses rj^ and 7700 — which are 

the ratios of and Kg amphtudes to a pair of pions, rjij = A{K^ 7r^7:^)/A{Kg vrV-') — so 
that 

2 1 

EK = -V+'^ + ^ Voo ■ (48) 

This is motivated by the so-called AI = ^ rule, i.e., that the amplitude to two pions in an / = 
state is approximately 20 times larger than the amplitude to two pions in / = 2. The definition 
of ek ensures that only the dominant hadronic amplitude contributes, so CP violation in decay 
gives negligible contribution to ek, and so ek ^ [i- — ^TnT{i=o)]/'^ (and 2Ree'ft' ?s 1 — \qK/PK\)- 
The phase of ek is independent of CKM elements (and of new physics), and is approximately 
7r/4 owing to to 2AmK « AFk- The current world average is \ek\ = (2.221 ± 0.008) x 10"^ [p^. 
Other observables measure \£k\ with less precision. Among these are the charge asymmetry in 
semileptonic decays, 6l, or decay-plane asymmetry in — > 7r"'"7r~e"'"e~ decays. 

Within the SM, ek is induced by AS = 2 transitions mediated by box diagrams. They are 
related to the CKM elements via a local hadronic matrix element, parameterized by the bag 
parameter Bk, 

{K'\{sj^{l - l,)df\K') = ^ mlflBKipi) , (49) 

where // is the renormalization scale, usually chosen to be 2 GeV. Note that Bk = 1 corresponds 
to the vacuum insertion approximation. The renormalization scale and scheme-independent bag 
parameter (sometimes denoted by Bk) is hereafter denoted by Bx, which differs from BxifJ-) 
by a series in that cancels the scale and scheme dependences (similar to B^ introduced in 



Sec. 3^). The kaon decay constant, fx = (159.8 it 1.5) MeV |45|, is extracted from the leptonic 
decay rate T{K~^ — > /x"'"f^). The most reliable prediction of Bx is obtained from LQCD. At 
present, calculations are performed assuming SU(3) symmetry and using mostly the quenched 
approximation. The world average is Bx = 0.79 ± 0.04 ± 0.09 [|^,||, where the first error 
combines statistical and accountable systematic uncertainties^ and the second is theoretical. 
Analytical approaches based on the large- A'c expansion of QCD find a smaller value for Bx in 
the chiral limit Large chiral corrections could be the origin of this difference. At present, 
Bx is the primary source of theoretical uncertainty in the SM prediction of Ex- 

Neglecting the real part of the off-diagonal element, M12, of the neutral-kaon mixing matrix. 



one obtains |21|. 



+ 2mtS{xc, xt) lm[V,sV:^VtsVt2 } , (50) 

where Amx = (3.490 ± 0.006) x IQ-^^MeV [||, and S{xi,Xj), with Xi = mf/m'^ {i = c,t), 
are Inami-Lim functions |p2[. The running quark masses are evaluated in the MS scheme. 



The term "accountable systematic uncertainties" denotes uncertainties that have been evaluated in a sys- 
tematic manner. Most of the experimental systematic uncertainties fall into this category. An example of an 
unaccountable (i.e., to a large extent arbitrary) theoretical uncertainty is the quenched approximation used in 
some LQCD calculations. 
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The parameters r/jj in Eq. (50) are next-to-leading-order QCD corrections to the Inami-Lim 



functions [IOC], where rycc is the parameter with the largest uncertainty. 



4.2 Direct CP Violation 

A nonzero value for the CP- violation parameter ej^, defined by 

1 



4 



{r]+- - mo) 



(51) 



establishes direct CP violation, because (X^o^o — X^+^-) /6. The experimentally convenient 
observable Re{e'j^/eK) is determined by the measurement of the double ratio 



r(i^o 



vrOvrO) 



ml 



TT ' TT 



Voo 



~ 1 — 6 Re 



(52) 



which cancels many experimental uncertainties. Evidence for direct CP violation was first seen 
by the NA31 Collaboration |lOll, and was firmly established by NA48 |0| and KTeV jlO^ , 
with the world average Re(e^/ei^) = (16.7 ± 1.6) x lO^''. 

Because Ke{e^/eK) is proportional to lm.{VtdVf*) = A^X^r] + ©(A^), it measures r]. However, 
the SM prediction of Re(e^/eft') suffers from large uncertainties, because the dominant gluonic 



and electromagnetic penguin contributions tend to cancel each other [104]. Detailed calculations 
at next-to-leading order ]105| , |106[ have been carried out for the coefficients of the two dominant 

(1/2) ('3/2') 

hadronic parameters, Bq (gluonic penguins) and (electroweak penguins), where the 

superscripts denote the isospin change in the K — > tttt transition. Because no reliable LQCD 
calculation for Bq^^'^^ is available, the measurement of Re(e^/eft:) is usually not used to constrain 
CKM elements. 



4.3 K+ 



TT'^uv and 



The E787 experiment at the Brookhaven National Laboratory (BNL) observed two events of the 
very rare decay 'k'^vV, resulting in B{K~^ — > tt^vu) = (l-57lQg2) ^ lO^"*^*^ [10'^]i which, 

owing to the small expected background rate (0.15 it 0.03 events), effectively excludes the null 
hypothesis. One additional event was observed near the upper kinematic limit by the successor 
experiment, BNL-E949 |10^ ]. They quote the combined branching fraction B{K~^ 7r~^i>V) = 
{lA7tlil) X 10-1°. 

In the SM, the branching fraction at next-to-leading order is given by ]ll09| 



■K VV] 



rK+- 



7^0e^ 



27r2 \Vus\'^sm^i 



'w 



xixt)VtdVZ + x.VcdV;, 



(53) 



where rj^+ = 0.901 corrects for isospin breaking ]110|. X{xt) = rjxXQ^Xt) is an Inami-Lim 
function for the dominant top-quark contribution, Xq, corrected by a perturbative QCD factor, 
rjx- The second term contains the charm-quark contribution. 



Expressing Eq. (53) in terms of the Wolfenstein parameters gives 



a L 



(54) 
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with fj = 1 + A2 + OiX^) and 

It provides an almost elliptic constraint in the (p, fj) plane, centered near (p = l,r/ = 0). The 
parameter k+ = A® [3a'^B{K~^ 7r'^e'^z^)]/[27r^ sin'* ^w) is such that B{K~^ tt^vu) depends on 
(j4A^)*, which is constrained by \Vcb\- The parameter Pq = -fc + ^Pc quantifies the charm-quark 
contribution, which until recently provided the largest theoretical error owing to the rric and 
scale dependences. The short-distance part, Pc = 0.38 it 0.04, is now calculated to next-to-next- 
to-leading order |111|. The long distance contributions of c and u loops are estimated to be 



5Pc = 0.04 lb 0.02 and are expected to be under good control |112]. 

In the SM, the decay K° — > TT^uu is theoretically clean, because the final state is almost 
completely CP-even |113||, so the decay proceeds dominantly through CP violation in the in- 



terference of decay with and without mixing |114, 115]. The amplitude is dominated by the 



top-quark contribution. The theoretical prediction of the branching fraction is given by | 109 



A5 



■X{xt) 



KLA''r^'X\xt) + 0{\^), (56) 



where kl = K+(ri^^ri^^)/(r^+r^+) = (2.12 ± 0.03) x 10"^° ]TT|, and rx^ = 0.944 accounts 
for isospin breaking [|110| . As for the charged mode, the branching fraction is proportional to 
(^A^)*, which would, however, cancel in the ratio of rates. The constraint in the (p, r?) plane 
from a future measurement of B{K^ tt^uu) would correspond to two horizontal bands at a 
certain value of ±|r/|. 

The interpretation of rare K decays involving lepton pairs is more complicated than for the 
similar B decays, because in most channels long-distance contributions are important. For 
example, — > Tr^£~^i~ is dominated by a virtual photon converting into £^£~ , whereas — >■ 
i~^£~ is dominated by virtual two-photon contribution. More interesting are the ir^i^i^ 
modes. Similar to ir^vV, the CP-violating amplitude dominates in — > vr'^e+e^, whereas 



there are comparable CP- violating and CP-conserving contributions in n^n^fi" [117]. A 

more precise measurement of B{Kg — > 7r^e~^e~) would allow the clean extraction of the relevant 
short-distance physics from K'^ — > TT^e~^e~. 

5 Unitarity Triangle Angle Measurements 

The observation of a nonzero Jarlskog parameter |14] in the SM leads to the requirement that 
the angles of the UT (and equivalently of all six unitarity triangles of the CKM matrix) be 



nonzero (modulo tt). The UT angles a, f3 and 7 (18) are all accessible from the B sector, albeit 
with different sensitivity and purity. Whereas the measurements of f3 (the leading experimental 
observable here is sin2/3) and 7, through B decays in charmonium and open charm, respectively, 
are theoretically clean, the measurement of a in charmless B decays relies on theoretical as- 
sumptions. Because the measurements of a and 7 involve interference with transitions governed 
by the small CKM matrix element V^h, they require larger data samples than when measuring 
sin2/3. 
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The experimental techniques to measure the UT angles also change radically from one to 
another. The measurements of a and /3 require mixing and therefore use neutral B 

mesons, whereas the measurements of 7 use interference between b ^ u and b ^ c decay 
amplitudes, and can be done with both neutral and charged B decays. 

In the following we neglect CP violation in B^ mixing, which has been searched for with both 
flavor-specific and inclusive B^ decays in samples where the initial flavor state is tagged. The 
current world average is \q/p\ = 1.0018 ± 0.0017 [[L18| , |11S(| , whereas the deviation from unity is 
expected to be \q/p\ — 1 ~ 0.0003 [120|, and around times smaller in B^ mixing. See Sec. |^ 



for an application of this search in the context of new physics studies. 
5.1 /3 from B Decays to Charmonium Final States 

In 6 — > CCS quark-level decays, the time-dependent CP violation parameters measured from the 
interference between decays with and without mixing are Sees = ~ilcP sin2/3 and Cccs = 0, to a 
very good approximation. The theoretically cleanest case if B ^ J/ipK^ where 

" ^ \y^J \y^J \y^J " ' ^''^ 

and so ImA^^o = S^^o = it sin 2/3 (see Sec. |1.4|) . The sign is from r]^j^o = =Fli the first 

S,L S,L S,L 

factor is the SM value of q/p in B^B^ mixing, the second is A/A, and the last one is px/qK- In 
the absence of K^K^ mixing there could be no interference between B^ — > ipK^ and B^ ipK^. 
BABAR pi and Belle pi have both used the rjcP = -1 modes J/i^K^, 'iIj{2S)K^, XciK° and 



rjcK^, as weh as J/^jKI, which has r]cP = +1, and J/il^ K*'^{892){^ K^'^'^), which is found from 
an angular analysis to have rjcP close to +1 (the CP-odd fraction amounts to 0.217 ±0.010 ||6^ ). 
In the latest result from Belle, only J/ipK^ and J/x/jK^ are used. The world average reads 

sin2/3 = 0.687 ± 0.032 , (58) 

giving for the angle (3 within [0,7r] the solutions (21.7l]^'2)° (68-3l|'3)°, where the first 
number is compatible with the result from the global CKM fit without the measurement of 
/?, {2AAtji)° and sin2/3cKM = 0.752 (cf. Sec. |). As expected in the SM, no direct CP 
violation has been observed in these modes. 

In b ccd quark-level decays, such as B^ J/iI)-k'^ or B^ D^*^D^*\ unknown contributions 
from (not CKM suppressed) penguin-type diagrams, carrying a different weak phase than the 
tree-level diagram, compromises the clean extraction of sin2/3. Consequently, they are not 
included in the sin2/3 average. 

5.1.1 Resolving the four- fold ambiguity 

Despite the agreement of Eq. (|5^ ) with the SM, it is still possible that, because of contributions 
from new physics, the correct value of /? is one of the three solutions not compatible with the 
SM. The measurement of the sign of cos2/? eliminates two of the solutions.^ P-meson decays 
to the vector- vector final state J/ipK*^, where three helicity states of the vector mesons mix 



The invariance /3 tt + /3 remains. It cannot be lifted without theoretical input on a strong phase 123|. 
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CP-even and CP-odd amplitudes, are also mediated by the b — > ccs transition. When the final 
state K*^ KgTT^ is used, a time-dependent transversity analysis can be performed allowing 



sensitivity to both sin2/3 and cos2/3 [124|. Here, cos2/? enters as a factor in the interference 
between the CP-even and CP-odd amplitudes. In principle, the sign of cos2/3 in this analysis 
is ambiguous owing to an incomplete determination of the strong phases occurring in the three 



transversity amplitudes. BABAR resolves this ambiguity by inserting the known variation [125| of 



the rapidly changing P-wave phase relative to the slowly changing 5-wave phase with the invari- 
ant mass of the Ktt system in the vicinity of the K*^{892) resonance. The result from BABAR, 
cos2/9 = 3.32 ILq'^q ± 0.27, although seemingly significant, exhibits a strongly non-Gaussian be- 
havior, and the confidence level for a positive value (as expected in the SM) is only 86% |127]. 



The result from Belle is compatible with zero |128| 



Bondar et al. [129| suggested another idea to resolve the four-fold ambiguity of /?. For quark- 
level transitions b — > cud, where the final state is a CP eigenstate (such as Dcpir^), the time 
dependence is given by the usual formulas, with the sine coefficient sensitive to sin2/3. When 
a multibody D decay, such as D Kg7:~^7r~ is used instead, a time-dependent analysis of the 
Dalitz plot of the neutral D decay allows the extraction of both sin2/3 and cos2/3. The Belle 
Collaboration [13C| has perfomed such an analysis, giving /3 = (16 it 21 it 12)°, which rules out 



the non-SM solution at 97% C.L.. Although HFAG has not attempted a combination of the two 
results on cos2/?, it is not premature to conclude that the SM solution is largely favored by the 
data. 

5.2 a. from Charmless B Decays 

Unlike P" J/'ipK^, for which amplitudes with weak phases different from the dominant tree 
phase are doubly CKM suppressed, multiple weak phases must be considered in most of the 
analyses of B decays to final states without charm. This complication makes the extraction 
of the CKM couplings from the experimental observables considerably more difficult, although 
richer. The decays most sensitive to a are B^ — > 7r"'"7r~, /O^vr^, and p~ (a measures vr — /5 — 7 
in the SM and in any model with unitary 3x3 CKM matrix) . The extraction of a in the presence 



of unknown penguin amplitudes requires an isospin analysis [131 1 for vrvr, pp, and a Dalitz-plot 
analysis [ |132| for p^vr^, which are briefly introduced below. Relying on flavor symmetries, in 
particular SU(2), does not represent a severe theoretical limitation. However, it certainly creates 
model-dependent uncertainties from flavor-symmetry breaking so that — neglecting statistical 
considerations — the measurement of a is not of the same quality as the measurements of sin2/3 
and 7. 

In naive factorization, there is a hierarchy between penguins in modes where the B^ decays to 
two pseudoscalar particles (such as vr+vr^), to a pseudoscalar and a vector particle (p^vr^), and 
to two vector particles (p+p^) [|133|| . This is due to the Dirac structure of the {V — A) x {V + A) 
penguin operators, which do not contribute when the meson that does not receive the spectator 
quark (the "upper" meson) is a vector, as in P" — > p~^7r~ and B^ — > p~^p~. Similarly, these 

^"The result for the strong phase is in agreement with the prediction obtained from s-quark helicity conserva- 
tion ITii. 
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Figure 5: Tree (left) and penguin (right) diagrams for the decay — > vr+vr 



operators interfere constructively (destructively) with {V — A) x (V — A) penguin operators 
when the upper meson is a pseudoscalar and the lower meson is a pseudoscalar (vector), as in 
vr^vr^ {B^ — > p~7r~^). This qualitative picture^^ reproduces quite well the experimental 

h0.41 
-O.IOj 



results. Using the value of a from the global CKM fit, the penguin-to-tree ratios are 0.23^^'^^ 



0.05+[];[!^, 0.03l[j;g^ and 0.10j:5];[!^, for vr+vr-, p+p-, p+ir- and p-7r+, respectively H, where 
penguin and tree correspond to the and amplitudes defined in Eq. (^) below. 

5.2.1 Relations Between Transition Amplitudes 

The general form of the B^ — > ■k~^tt~ decay amplitude, accounting for the tree and penguin 
diagrams that correspond to the three up-type quark flavors (u, c, t) occurring in the W loop 
(see Fig. |5|), reads 

A+- = AiB"^ ^ vr+TT-) = VudV:i,Mu + V^V:,M, + VtdVtlMt , (59) 

and the CP-conjugated amplitude can be similarly written. Using unitarity, this can be written 
asi2 

A+- = VudV:,T+- + VtdV;,P+- , (60) 

with r+- = Mu-Mc and P+~ = Mf- M^. We refer to the T+- and P+~ amplitudes as "tree" 
and "penguin", respectively, although it is implicitly understood that both of them receive 
various contributions of distinct topologies, which are mixed under hadronic rescattering. The 
time-dependent CP-violating asymmetry is given by Eqs. (24) and (^), with the CP-violation 



parameter A^+^- = {q/p)[A'^ /A^ ), since vr+vr" is a CP eigenstate with eigenvalue +1. In 
the absence of penguin contributions A^+^- = e^*", and hence 5*^+7^- = sin2a and C^+^- = 0. 
In general, however, the phase of X-j^+t^- is modified by the interference between the penguin and 
the tree amplitudes, which can also lead to direct CP violation. It is customary to define an 
effective angle agg that incorporates the penguin-induced phase shift A^+^- = |A^+^-| e^*"<=ff, so 
^^Because iV ~ A) x {V + A) operators are power suppressed in the QCD factorization approach, the naive 



hierarchy may receive large corrections [134| 

12 



In Eq. (pCj), the c-quark loop is eliminated using the unitarity relation. In Charles et a/.'s definition, this 
corresponds to the c-convention. Note that the particular choice of which amplitude to remove in the definition 
of a total transition amplitude is arbitrary and does not have observable physical implications. 
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that the sine coefficient becomes S^+^- = y 1 — C^+^- sin2aeff. Twice the effective angle Oeff 
corresponds to the relative phase between e~'^'^^A'^ and A'^ . 

The isospin invariance of the strong interaction relates the amplitudes of the various B vrvr 
decays to each other. A pair of pions in a i? — > vrvr decay must be in a zero angular momentum 
state, and then, because of Bose statistics, they must have even isospin. Consequently, tt^tt^ is 
in a pure isospin-2 state, free of QCD penguin amplitudes, which contribute only to the isospin- 
final states. Gronau & London showed that the measurements of rates and CP-violating 
asymmetries of the charged and two neutral tttt final states provide sufficient information to 
extract the angle a as well as the various tree and penguin amplitudes yi31| . Unfortunately, 
this isospin analysis is plagued by a 16-fold ambiguity for a S [0,27r], i.e., 16 (not necessarily) 
different values of a reproduce the same set of observables. 

Using the same convention as in Eq. (|60|), one can write for the remaining B tttt decays 



^/2A00 ^ V2AiB' ^ vrVO) = V^T^' + UtdU;,P°° , (61) 

and the CP-conjugated modes can be similarly written. The subscript C stands for a color- 
suppressed amplitude, and the subscript EW stands for the electroweak penguin amplitude 
contributing to vr^vr''. The Pew contributions are usually neglected in the amplitude fits per- 
formed by the experiments, so that the A/ = 3/2 decay P"*" 7:~^tt^ is mediated by a pure tree 
amplitude. Applying the isospin relations [|131[| 



A+^ = ^A+- + A'\ and A+' = ^A+- + A^^^ (62) 



the amplitudes (61) can be rearranged and, e.g., T^^ and P^^ eliminated. Isospin-breaking 
effects have been examined in Refs. [138| and are expected to be smaller than the current 
experimental precision on a. 

Applying the isospin relations reduces the number of unknowns in the B — > vrvr isospin analysis 
to six, which aligns with the number of experimental observables: three branching ratios, one 
mixing-induced CP asymmetry (5^+^-), and two direct CP asymmetries (C^r+Tr- and C^o^o). 
Owing to the missing tracks in most of the P" vr'^vr" events, 5*^07^0 cannot be measured in 
the foreseeable future. In absence of isospin breaking, direct CP violation cannot occur in 
B~^ — > vr"''vr''. A small value for the branching fraction to vr^vr*^ would indicate that the penguin 
contribution cannot be large |139| . This observation led to the development of isospin bounds 



on Aa = a — agg |140, 141]. Applying the optimal bound is equivalent to constraining a using 



the isospin relations without using the measurement of C^o^o . 
5.2.2 Isospin Analysis of S — > tttt Decays 

All six observables used in the isospin analysis have been measured in P — > vrvr decays. The 
past disagreement between BABAR and Belle for the CP asymmetries in B^ — > vr^vr" (Belle 



^''The bulk part of the shift of a due to Pew can be estimated model independently |135 
to -2.1° H. 



136 



137 , and amounts 



CP violation and the CKM matrix 



34 




Figure 6: Left: confidence level for the UT angle a obtained from the B — > vrvr isospin analysis 
under various conditions: full isospin analysis (dark shaded), isospin analysis without using 
C^o^o (light shaded), and full isospin analysis assuming perfect knowledge of the parameters 
Sj^+T^- and C^+TT- (dashed line). Right: confidence level for the UT angle a from the isospin 
analyses of B — > tttt (dotted) and B pp (dashed- dotted), and the Dalitzplot analysis of B pir 
(dashed). The shaded region shows the combined constraint from all three measurements jji]/. 



saw evidence for both mixing-induced and direct CP asymmetries, whereas BABAR did not) has 
been alleviated recently. The world averages for the CP asymmetries and branching fractions 
measured by BABAR, Belle and CLEO (for the branching fractions) have been compiled by 
the HFAG [142, 66 1. A large B^^/B^~ ratio is found, indicating significant penguin or color- 
suppressed tree contributions, and representing a challenge for the theoretical understanding of 
the decay dynamics. This ratio implies only a loose upper bound |Aq| < 40° at 95% C.L., if 
C^o^o is not used [51|. 

The confidence level for a from B vrvr is shown in the left plot of Fig. |6[ Three cases 
are distinguished: a) the full isospin analysis using the present world averages (shaded), b) the 
incomplete isospin analysis without C^o^o (solid line), and c) the full isospin analysis assuming 
perfect knowledge of S^+tt- and C^^+j^- (dashed line). The present constraint on a from tt^tt" 
is weak, and future improvement will be driven almost entirely by the measurement of C^o^o. 



5.2.3 Isospin Analysis of B — ^ pp Decays 

Similar to B ^ tttt, the angle a can be extracted from the measurement of time-dependent CP 
asymmetries in B^ — > p^p^ , together with an isospin analysis including the full B ^ pp system, 
to determine the penguin pollution. The specific interest of this mode is the potentially small 
penguin contribution, which is suggested by the smallness of the experimental upper bound on 
the B^ — > p^p^ branching fraction when compared with the branching fractions of B^ P~^P~ 
and B'^ P^P^ [143, 66|. In addition, both direct and mixing-induced CP-violating asymmetries 
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of the — > p^p^ decays are experimentally accessible once the decay has been seen. 

The analysis is complicated by the presence of three helicity states for the two vector mesons. 
One corresponds to longitudinal polarization (^o) and is CP-even. Two helicity states {A±) are 
transversely polarized and are admixtures of CP-even and CP-odd amplitudes. Because helicity 
flips are suppressed in these decays, the relative proportions of the helicity amplitudes in naive 
factorization are expected to be |144] 

.2 



rub \ mb J 



(63) 



The dominance of the longitudinal polarization has been confirmed by experiment in the mea- 
sured B ^ pp decays. The isospin analysis can therefore, without significant loss of sensitivity, 
be applied to the longitudinally polarized states only.^^ 

The confidence level obtained from the B ^ pp system for a is given in the right plot of 
Fig. ^ (dashed-dotted line). The reduction of B{B^ — > p^p^) by the latest BABAR measurement 
improved the consistency of the isospin analysis, which now exhibits a closed isospin triangle, 
and deteriorated somewhat the resulting constraint on a. Choosing the solution consistent with 
the SM gives app = (100 ^20)° lUII- Also shown in the plot is the prediction from the CKM 
fit, where the direct measurement of a has been excluded from the fit. It exhibits excellent 
agreement with the measurement. The 95% C.L. bound on Aa is |Aa| < 17°, reflecting a small 
penguin pollution. 

5.2.4 Dalitz Plot Analysis of B — > pvr Decays 

Unlike vr^vr" and p^p^, the final state p^ir^ is not a CP eigenstate, and four flavor-charge 



configurations {B^{B ) — > /> vr^) must be considered. The corresponding isospin analysis |14(;] 



is unfruitful with the present statistics, because two pentagonal amplitude relations with a 



total of 12 unknowns have to be solved. However, as Snyder and Quinn [132| pointed out, 
the necessary degrees of freedom to constrain a without ambiguity (except for the irreducible 
a — > a + vr), can be obtained by including in the analysis the (known) strong-phase modulations 
of the interfering p resonances in the Dalitz plot. 

The BABAR |147] ] Collaboration has performed such an analysis, where the interfering B^ 
decay amplitudes into p'^ir^, p^ir^ and p^n^ have been included. To avoid mirror solutions for 
the fit parameters and non-Gaussian observables, BABAR fits 16 interdependent coefficients that 
are bilinears of the charged and neutral p form factors. Each coefficient is related in a unique 
way to the quantities of interest, the tree-level and penguin-type amplitudes, and a. These 
are obtained in a least-squares fit to the measured coefficients. Systematic uncertainties are 
dominated by the lack of knowledge of the tails of the p form factor, where higher p excitations 
occur and precise B decay data are scarce. Ignoring the mirror solution at a + 180°, BABAR 



Owing to the finite width of the p meson, 1 = 1 contributions can occur in B ^ pp decays [f45]. Although no 
exphcit calculation is available, one expects these effects to be of order {Vp/nipY ~ 4%. At present the experiments 
also neglect effects from the interference with the radial excitations of the p, with other tttt resonances, and with 
a nonresonant component. The latter two effects are included in the systematic errors quoted by the experiments. 

^^The full set of bilinears comprises 27 observables out of which BABAR neglects nine related to CP asymmetries 
in P^'tt'', which has only a small event yield. 
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finds OpTT = (113 =b 6)°, whereas only a marginal constraint on a is obtained beyond two 
standard deviations. 

The confidence level for a obtained from the B ^ pn Dalitz analysis is given in the right 
plot of Fig. ^ (dashed line). Also shown is the combined constraint from all three measurements 
(shaded), giving 

a={l00tfj (64) 

for the SM solution. This measurement is in agreement with the expectation acKM = (97l5g)° 
from the global CKM fit (see Table | on p. ||) |5|. 



5.3 7 from B Decays to Open Charm 

The golden methods to determine 7 at the i?-factories utilize the measurement of direct CP 
violation in DK^ decays, where the neutral D meson can be both D° and (and 

where also stands for D*^). The corresponds to the leading b transition, whereas the 
is produced by a CKM- and color-suppressed b ^ u transition. If the final state is chosen so 
that both L»° and can contribute, the two amplitudes interfere, and the resulting observables 
are sensitive to the UT angle 7, the relative weak phase between the two B decay amplitudes. 

Among the many methods that exploit this interference, the experiments concentrate on the 
reconstruction of the neutral D in a CP eigenstate (GLW) [ |148| |, in other final states common to 
D'^ and such as K^tt^ (ADS) |149f| , or in the self-conjugate three-body final state K^tt^tt^ 
(GGSZ) [ p.50| . All variations are sensitive to the same B decay parameters and can therefore 
be treated in a combined fit to extract 7. For the GLW analysis, the relations between the 
experimental observables, the rates Rcp± and charge asymmetries A(jp± of the CP-even and 
CP-odd modes, and the model parameters are 

2 r A 2r B sm.5 B sm.'y 

Rcp± = l + rB±2rBCosdBCOs-j, Acp± = ± , (65) 

ncp± 

where 6b is the (CP-conserving) relative strong phase between the interfering amplitudes, and 
r^j*^ = \A{B+ D^(*^K+)/A{B+ D°'^*^K+)\ is the ratio of the B decay amplitudes. The 
same relations apply for the GGSZ analysis, where the variations of additional strong phases 
and amplitude ratios throughout the Dalitz plot of the three-body decay are extracted from high 
statistics D decay data and modeled by the experiments. For the ADS analysis two additional 
unknowns (phase and amplitude ratio) appear because the neutral D mesons do not decay 
to CP eigenstates. High-statistics charm decays can be used as external constraints for the 
amplitude ratio. Moreover, only two observables (rate and asymmetry) per B decay mode exist 
for ADS. The above equations show that the feasibility of the 7 measurement depends crucially 
on the size of (expected to be roughly ~ 0.1, if naive color suppression holds). Instead 
of fitting the model parameters (7, r^, 6b) directly, the experiments proceed by determining 
the Cartesian coordinates {x±,y±) = {tb cos{6b ± 7), sin(5B it 7)), which have the advantage 
of being Gaussian distributed, uncorrelated and unbiased (the parameter is positive definite 
and hence exhibits a fit bias toward larger values, if its central value is in the vicinity of zero 
within errors), and simplify the averaging of the various measurements. 
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Results are available from both BABAR and Belle on the GLW analysis using the decay modes 



B+ DK+, B+ D*K+ and S+ DK*+ pi. For ADS pi only BABAR uses all three 



decays. Here, the decays D* — > Dtt^ and D* have a relative strong-phase shift and 



must be studied separately |153|. The world averages for the GLW and ADS measurements 
are provided by the HFAG . So far none of these measurements have obtained a significant 
signal for the suppressed amplitude. 

Both experiments have studied the GGSZ analysis using the same B decay modes as for the 
GLW and ADS analyses, and using in all cases the subsequent decay D^{D^) K^'k^-k~ . Both 
experiments perform a full frequentist statistical analysis to determine the physical parameters. 
Exploiting the three-body Dalitz plot enhances the sensitivity to the suppressed amplitude 
compared to the (quasi-)two-body decays used for GLW and ADS. This leads to a significant 
determination of tb and hence 7 with this method. Using the same three modes 
the experiments report the GGSZ measurements 7 = (67 ± 28 ± 13 ± 11)° (BABAR ||T5l ) and 
7 = (53l]^|ib3ib9)° (Belle [|155| ]), where the first errors are statistical, the second are systematic, 
and the third are due to the Dalitz plot model. 

The CKMfitter Group has performed a combined analysis of all GLW, ADS and GGSZ results 
using a frequentist statistical framework. Including the former two methods leads to a reduction 
in the average rs values of all modes, which are rB^DK) = 0.10 it 0.03, rB{D*K) = 0.10 lg o| 
and rB{DK*) = 0.111051 [^l|. All values are in agreement with the naive expectation from 
CKM and color suppression. The smaller tb values deteriorate the constraint on 7 obtained 



from the GGSZ analyses. The preliminary overall 7 average is [51| 

7 =(6211)°, (66) 



which is in agreement with the prediction 7ckm = (6011)° from the global CKM fit (where 
the direct 7 measurement has been excluded from the fit) Q|. The GLW, ADS and GGSZ 
confidence levels for 7 as well as the combined constraint are shown in Fig. 0. Also shown is the 
expectation from the global CKM fit. 



5.4 2/3 + 7 from B Decays to Open Charm 

Similar to the decay B^ — > p^vr^, which is not a CP eigenstate but sensitive to a because 
both final states can be reached by both neutral B flavors, interference between decays with 
and without mixing can occur in B^ D^*^^Tr^{p^). A time-dependent analysis of these 
decays is sensitive to sin(2/3 + 7), because the CKM-favored 6 — > c decay amplitude interferes 
with the CKM-suppressed b ^ u decay amplitude with a relative weak-phase shift 7. In these 
b — > c{ud), u{cd) quark-level transitions no penguin contributions are possible, because all 
quarks in the final state are different. Hence there is no direct CP violation. Furthermore, 
owing to the smallness of the ratio of the magnitudes of the CKM-suppressed and CKM-favored 
amplitudes (denoted rj), to a very good approximation, C/ = —Cj = 1 (using 
/ = D^*^^h^ h = Tr,p), and the coefficients of the sine terms are given by 



Sf = -2rf sin(2/? + -f-6f), Sj= -2r/ sin(2/3 + 7 + dj), 



(67) 
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Figure 7: Left: confidence level for the UT angle 7 from the measurement of direct CP asym- 
metries in B ^ decays (left). Shown are the individual results for the combined GLW 
and ADS analysis (dashed- dotted), for GGSZ (dashed), and the combination of both analyses 
(shaded). Right: adding to these the measurement of 2(3 + ^ from time- dependent analyses of 
neutral B decays to open charm, together with the precise sin2/3 measurement. A theoretical 
error of 15% has been used to account for SU(3) breaking. The plots are taken from Ref. J^/ . 



where (5/ is the strong-phase shift due to final-state interaction between the decaying mesons. 
Hence, the weak phase 2/3-1-7 can be cleanly obtained from measurements of Sj and Sj, although 
external information on at least rj or (5/ is necessary. It can be obtained from the corresponding 
flavor-tagged branching fractions, or from similar modes that are easier to measure. These can be 
ratios of branching fractions of the charged B^ Z)(*)+7r^ to the neutral CKM-favored decay, or 
ratios involving self-tagging decays with strangeness like B^ Di*^^7r-. Corrections for SU(3) 
breaking in the latter case generate a significant theoretical uncertainty, which is in general hard 
to quantify. Naively, one can estimate r^*) ~ | V^'^Kib/Kid^cbl — 0.02. At present, the most 
precise estimate of obtained from SU(3)-corrected ratios are |156[| r^i^ = 0.013 ± 0.002 and 
rD*n = 0.019 lb 0.0029, where a 15-30% theoretical uncertainty is usually added to account for 
SU(3) breaking. 

The constraint on 2/3 -|- 7 obtained using these amplitude ratios (assuming a 15% uncertainty 
on SU(3) breaking), and the world averages of the experimental observables |157| , |66[ , is depicted 
by the dotted line in the right plot of Fig. 0. The measured values for sin2/3 (and cos2/3) have 
been used to infer 7 from this. The shaded curve gives the overall constraint on 7 obtained 
when averaging these results with those given in Sec. 5^. The 2/3-1-7 rneasurements help to 
exclude large values of 7. 



5.5 Weak-Phase Information from Direct CP Violation in B Decays 

The KM mechanism causes "direct" CP violation in the decay, as soon as at least two amplitudes 
with different strong and weak phases contribute. Because virtual loops are present in all meson 
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decays, "some" (possibly unobservable) amount of direct CP violation always occurs. Owing 
to the large weak phases arising in B decays, direct CP violation should be more prominent 
here than, e.g., in the kaon system. This has been confirmed by the measurement of the direct 
CP asymmetry Ak+^- = -0.115 ± 0.018 in B° K+tt' decays Evidence for direct CP 

violation in neutral B decays also exists for B'^ tt^tt^ (S.Tcj significance) and B^ — > p^vr^ 
(3.4(T significance) p6|. Recently, the first evidence for direct CP violation in charged B decays 



emerged from the mode B~^ — > with a charge asymmetry of 0.31 IgJ^ [f59|. With the 



-0.11 

data samples at the i3-factories increasing, we expect the discovery of more and more rare-decay 
modes with significant CP violation in the decay. 

From the point of view of the weak-phase extraction, the required conspiracy between compet- 
ing amplitudes of similar size and the occurrence of strong phases, represent serious obstacles. 
A reliable and model- independent calculation of direct CP violation is not possible at present, 
and estimates based on factorization are plagued by large uncertainties. However, flavor sym- 
metries in particular isospin can be exploited to (essentially) assess model independently direct 
CP violation. In B decays to vrvr, pvr and p/3, the measurements of direct CP-violating asymme- 
tries (independent of whether they are compatible with zero or not) are essential inputs to the 
isospin analyses. In the Kt: system the corresponding isospin analysis used to extract 7 |146| is 



fruitless at present, and affected by possibly large isospin-breaking corrections from electroweak 
penguins, which cannot be taken into account model independently as is the case in the vrvr and 
pp isospin analyses. 

Although a quantitative prediction is difficult, direct CP violation can be a powerful probe 
for new physics in decays where negligible CP asymmetries are expected. This is the case for all 
B decays dominated by a single decay amplitude. Prominent examples are penguin-dominated 
decays, such as 6 — > 57 or B K^*k+e-, where a significant nonzero direct CP violation would 
unambiguously indicate new physics. 

6 The Global CKM Fit 

The global CKM fit consists of maximizing a likelihood built upon relevant experimental mea- 
surements and their SM predictions, which depend on the parameters of the theory. Some of 
these parameters, such as quark masses or matrix elements, are experimentally or theoretically 
constrained, whereas others are unknown. These unknowns contain the four Wolfenstein pa- 
rameters, but also, for instance, hadronic quantities that occur in the determination of the UT 
angles a and 7. To avoid uncontrolled biases, the unknown parameters should be treated as free 
parameters of the theory in a frequentist statistical sense. Such an approach has been adopted 
by the CKMfitter Group [|l6^,|ll, whose results we use here. 



Many of the observables presently used to constrain the CKM matrix depend on hadronic 
matrix elements and hence are subject to significant theoretical uncertainties. The treatment of 
these theoretical uncertainties has a large amount of arbitrariness, and different choices lead to 



different results. The CKMfitter Group advocates the Rfit approach [160,13]. In this approach 
the ranges spanned by the model-dependent theoretical errors are scanned to maximize the 
agreement between theory and experiment, which corresponds to a minimization of the exclusion 
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confidence level at a given point in parameter space that is to be tested. Other approaches exist 



in the literature |161, 162]. 



6.1 Fit Inputs and Initial Tests of Unitarity 

What is often termed the "standard global CKM fit" includes only those observables for which 
the SM predictions can be considered quantitatively under control, and that lead to a significant 
constraint on the CKM parameters. These comprise the measurements of \Vud\, {Vusl, \Vub\, 
iKfel) I^kI, Amrf, Anis,^^ sin2/3, cos2/3, a, 7 and B{B~^ — > t'^Ut). The inputs used by the fit are 
given m Table |l[ 

Using the independently measured moduli of the CKM elements mentioned in the previous 
sections, the unitarity of the CKM matrix can be checked. We obtain |14dp + |Vusp + |Kifep — 1 = 
-0.0008 ± 0.0011 (first row), \Vcd\'^ + \Vcs\'^ + \Vcb\^ - 1 = -0.03 ± 0.18 (second row), and 
l^udP + {Vcdl"^ + |Vtdp — 1 = 0.001 lb 0.005 (first column). The sum in the second column, 
iKtsP + + l^sP is practically identical to that in the second row, as the errors in both 

cases are dominated by \Vcs\- For the second row, a more stringent test is obtained from the 
measurement of Sncdsfel^jP minus the sum in the first row, giving \ Vcd\'^ + \ Vcs\'^ + iVctp = 
1.000 ± 0.026. 

It is an unexpected success of the i?-factory experiments that, beyond the precision mea- 
surement of sin2/3, measurements of a and 7 became available. Remarkably, the most powerful 
determinations of the latter two angles involve measurements that were not considered before 
the data became available. Choosing the SM solution for each of the angles, their sum reads 

a + /3 + 7= (l86li?)° , (68) 

which corresponds to a closed triangle verifying the three-generation KM mechanism. 

Various theory parameters enter the SM predictions that relate the measurements to the 
fundamental CKM parameters determined by the fit. Those with the largest uncertainties are 
hadronic matrix elements obtained with LQCD, quark masses determined using LQCD and 
QCD sum-rule analyses, and coefficients calculated in perturbative QCD. The values used for 
these parameters in the CKM fit are also given in Table ^ 



6.2 Fit Results 

The global CKM fit results in the Wolfenstein parameters are given in Table (top four lines). 



The goodness-of-fit, evaluated with the use of Monte Carlo methods |13|, is satisfying, so one 
can assume the correctness of the SM and proceed with the metrology of the CKM parameters 
and the related quantities. It is interesting to compare the direct measurements of the UT 
angles with the predictions from the CKM fit that do not include these measurements. The 
measurement of sin2/3 outperforms by far the indirect prediction. Also, the precision of the 
a measurement exceeds that of the CKM prediction so that its inclusion in the fit improves 
the knowledge of the apex of the UT. For 7, the experimental precision is not yet sufficient to 
The results shown in this section include the new D0 and CDF analyses of B°B° mixing ||,|6). 
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Errors 



Parameter 


Value ± Error(s) Reference 


GS 


TH 


\Vud\ (/3-decays) 
\Vus\ (Kiz and Kf,2) 
\Vub\ (inclusive) 
\Vub\ (exclusive) 
\Vcb\ (inclusive) 
\Vcb\ (exclusive) 


0.97377 ± 0.00027 
0.2257 ±0.0021 
(4.45 ± 0.23 ± 0.39) X 10-3 | 

(3.94 ±0.28 ±0.51) x lO'^ 
(41.70 ±0.70) X 10-3 
(41.2 ± 1.7) X 10-3 


|98| 

il 

11 

il 
il 


i 




★ 
* 


\eK\ 
Arud 

AlTls 


(2.221 ± 0.008) X 10-3 1 

(0.507 ±0.004) ps-i 

Amplitude spectrum (incl. D0), and 
CDF log£ parameterization 


n 


* 




LiiiiiJ^ 1 CCI 


0.687 ± 0.032 


il 








B +10 -/o 
BO ^ (p7r)0 ^ 3^ 


Inputs to isospin analysis 
Inputs to isospin analysis 
Time-dependent Dalitz plot analysis 


il 

163 


★ - 

* - 


B- Di.*) k'^*)- 
B- d'-*'^K'^*^~ 

B- 


Inputs to GLW analysis 
Inputs to ADS analysis 
GGSZ Dalitz plot analysis 


il 
il 
il 








B{B- T-Vr) 


Experimental likelihoods 


il 








mc{mc) 
fnt{mt) 

fK 


(1.24 ± 0.037 ± 0.095) GeV 
(162.3 ± 2.2) GeV 
(159.8 ± 1.5) MeV 


il 

164 

n 








Bk 

Vet 

Vti 

r?B(MS) 
Bd 


0.79 ±0.04 ±0.09 
0.1176 ±0.0020 

0.47 ±0.04 
0.5765 ±0.0065 
0.551 ±0.007 
(191 ± 27) MeV 
1.37±0.14 
1.24 ±0.04 ±0.06 


il 
il 

100 
100 

n 

il 






★ 
★ 
* 



anticorrelated theory error with /s^ \/ Bd 



Table 1: Inputs to the standard CKM fit. If not stated otherwise, when two errors given, the 
first is statistical and accountable systematic and the second stands for theoretical uncertainties. 
The last two columns indicate whether errors are treated Gaussian (GS) or theoretical (TH) 
within Rfit JT^/ . Table taken from Ref. j57|/. 

improve the CKM metrology. It is remarkable that the prediction of sin2/3s has an absolute 
precision of 0.002 (owing to its CKM suppression), so a measurement of it at the LHC will 
represent a very sensitive probe for physics beyond the SM. Also given are the sides of the UT, 
and the angles of the standard parameterization of the CKM matrix. 

Once the Wolfenstein parameters are fit, determining fully consistent confidence levels for all 
related observables is straightforward. It is also possible to derive constraints on the matrix 
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Observable 


central it C.L. = \(7 


-1- P T — 9^ 


It = 00 


A 
A 

P 
fj 


PI 99 79 +0.0010 

0.809 iroll 

o.i97iroro 

n QQO +0.019 


+0.0020 
-0.0020 
+0.029 
-0.028 
+0.050 
-0.087 
+0.047 
-0.037 


+0.0030 
-0.0030 
+0.044 
-0.042 
+0.074 
-0.133 
+0.075 
-0.057 


J [10-5] 


o nr:+0.18 
■J-^^ -0.18 


+0.45 
-0.36 


+0.69 
-0.54 


sin2a 

sin2Q (meas. not in fit) 
sin2/3 

sin2/3 (meas. not in fit) 
a (deg) 

a (deg) (meas. not in fit) 
a (deg) (direct meas. only) 
13 (deg) 

P (deg) (meas. not in fit) 
P (deg) (direct meas. only) 
7 ~ (deg) 

7 ~ 5 (deg) (meas. not in fit) 
7 ~ 5 (deg) (direct meas. only) 
/3s (deg) 
sin(2/3,) 


"•^^ -0.15 

-0.23t°:f6 

n 71 f; +0.024 
^•'-■-D -0.024 

7^9+0.057 
\J.iOZ, _o.o35 

97.311:5 
96.51^6% 
100.2 +I5g0 

22.86 ll^o 
24.4 tij 
21.7tiii 
59.8 tt? 

59.8 t^l 
62 1| 

1 Q^tr +0.061 
l.U40_o,o57 

U.UODO _o.0020 


+0.49 
-0.28 
+0.72 
-0.32 
+0.048 
-0.049 
+0.105 
-0.073 
+8.7 
-14.0 
+9.9 
-21.2 
+22.7 
-18.2 
+2.03 
-1.97 
+5.1 
-3.0 
+2.6 
-2.4 
+13.9 
-7.9 
+14.1 
-8.0 
+61 
-39 
+0.151 
-0.114 
+0.0053 
-0.0040 


+0.71 
-0.42 
+0.83 
-0.45 
+0.074 
-0.075 
+0.135 
-0.112 
+13.7 
-20.7 
+14.6 
-25.3 
+32.0 
-28.1 
+3.22 
-2.93 
+6.9 
-4.5 
+4.1 
-3.6 
+20.8 
-12.1 
+21.0 
-12.3 
+100 
-54 
+0.238 
-0.177 
+0.0083 
-0.0062 


Ru 
Rt 


0.391 iroii 

n 070 +0.033 
U.O^Z _o 028 


+0.031 
-0.029 
+0.095 
-0.054 


+0.049 
-0.044 
+0.143 
-0.082 


Am^ (ps"!) (meas. not in fit) 
ArUs (ps""*^) 

Am^ (ps"l) (meas. not in fit) 
Ek [IQ-'^] (meas. not in fit) 


0.394 Ig-.Sr^ 

17.34 Ig-.fo 
21.7t5;9 

2.46l°o:i 


+0.219 
-0.132 
+0.65 
-0.35 
+9.7 
-6.8 
+1.05 
-1.05 


+0.361 

-0.162 

+0.78 

-0.49 

+13.1 

-9.1 

+1.50 

-1.20 


Jb^ (MeV) (lattice value not in fit) 
^ (lattice value not in fit) 
Bx (lattice value not in fit) 


18311° 
1 061 +0122 

J--tlDl _o.047 
n 799+0.251 
U. (ZZ _o.084 


+21 
-20 

-0.083 
+0.348 
-0.157 


+34 
-28 

1 n QO/i 
+U.OZ4 

-0.119 

+0.461 

-0.216 


rric (GeV/cP) (meas. not in fit) 
nit (GeV/c^) (meas. not in fit) 


U-ol-0.36 

1501^7 


+1.08 
-0.36 

+ / 

-35 


+1.23 
-0.81 

-46 


BiKl^TT^uu) [10-11] 
B{K+ TT+uV) [10-11] 


9 CO +0.48 

7 c;+1.8 
' -"3-2.0 


+1.01 
-0.68 
+2.5 
-2.4 


+1.53 
-0.93 
+3.2 
-2.7 


B{B+^T+iy^) [10-5] 

B{B+ T+u^) [10-5] (direct meas. only) 

B{B^^^i+v^) [10-7] 


n 7+1-3 
' -1.3 

10.411? 

A 09+0.58 

^■-J^ -0.57 


+2.9 
-2.5 
+3.0 
-5.1 
+1.27 
-1.12 


+4.6 
-3.6 
+3.0 
-7.1 
+2.05 
-1.62 



Table 2: Numerical results of the global CKM fit (I) /]5l|/. The errors correspond to one, two 
and three standard deviations, respectively. 
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Observable central ± C.L. = la ± C.L. = 2(7 ± C.L. = 3ct 



\v ^ 


1 
1 






07383+° "™^'* 


+0.00047 
-0.00047 


+0.00071 
-0.00071 


\V 
1 'us 


1 






r, 99 79 +0.0010 
U.^^(^„0.0010 


+0.0020 


+0.0030 


1 






-0.0020 


-0.0030 


\V K 


1 fio-^i 

1 L^'J J 






o 09 +0.15 
"^■"^-O.IS 


+0.31 
-0.29 


+0.49 
-0.44 


\V h 
\ * ub 






not in nf 1 


3 64+S J? 

>J.U4:_0.18 


+0.39 
-0.36 


+0.60 
-0.55 


\Vcd 


1 






r, 9971 9 +0.00099 
U.zz / IZ _o. 00103 


+0.00199 
-0.00205 


+0.00300 
-0.00307 










97297 +° °™?^ 

u.ai -0.00023 


+0.00048 


+0.00071 


1 






-0.00047 


-0.00071 


\Vcb 


1 [10-^] 






41 79+!!-?3 

( »_o.63 


+ 1.26 

-1.27 


+ 1.89 

-1.90 


\Vcb 


1 [10-^] 1 




not in fit 1 

XXW U XXX XX L / 


44 q+i ? 
^^■^ -2.8 


+2.4 
-5.7 


+3.8 
-7.7 


\Vtd 


1 [10-^] 






-0.29 


+0.92 
-0.57 


+ 1.38 
-0.86 


\Vts\ 


[10-3] 






41 13+!??o 

"±±.±0 _o.62 


+ 1.25 
-1.24 


+ 1.87 
-1.86 


\Vtb\ 








999119 

u.aja±± J -0.000027 


+0.000052 


+0.000078 








-0.000054 


-0.000082 


\VtdlVts\ 






2011+" ™^^ 
u.zuii -0.0065 


+0.0230 
-0.0127 


+0.0345 
-0.0195 


\Vud 


KTJ [10- 






' ^ -0.14 


+0.30 
-0.29 


+0.48 
-0.43 


arg 


[VudV:,] 


(deg) 




59 8+!-^ 


+ 13.9 
-7.8 


+20.9 
-12.1 


arg 




(deg) 




1 043+S-Sei 

i.vio —0.057 


+0.151 
— 0.114 


+0.238 
— 0.176 


\VcdV:,\ [10- 






9.49t°:}^ 


+0.30 
-0.30 


+0.45 
-0.45 


arg 


[-Vcdv:^] 


(deg) 


1 


ij.yjoo\) -0.0020 


+0.0050 
-0.0040 


+0.0077 
-0.0060 


\Vtd 


v:b\ [10- 






97 +0.33 
-0.29 


+0.93 
-0.57 


+ 1.38 
-0.85 


arg 


[VtdV;,] (deg) 




-22.84 t^.-oo 


+ 1.98 
-2.02 


+2.93 
-3.21 


sin 6 


'l2 






n 0070 +0.0010 
u.zzrz_o.ooiO 


+0.0020 
-0.0020 


+0.0030 
-0.0030 


sin 6 


'l3 [10-3 


] 




q 09+0.15 
"^■»^-0.15 


+0.31 
-0.30 


+0.49 
-0.44 


sin 6 


'23 [10-3 


] 




41.78 ^SJJ3 


+ 1.26 

-1.26 


+ 1.90 

-1.89 



Table 3: Numerical results of the global CKM fit (II) |5l|/ . The errors correspond to one, two 



and three standard deviations, respectively. 



elements of neutral K and B mixing, which in turn are in agreement with the LQCD results. 
The determination of through the Am^ measurement in the CKM fit is not yet competitive 
with the LQCD determination, owing to the currently insufficient precision of the 7 and a 
measurements. It represents a remarkable success of flavor physics to correctly reproduce the 
scales of the charm and top masses through FCNC loop processes only. Also listed in Table || are 
the predictions for the rare kaon decays to pions and neutrino-antineutrino pairs (including the 
next-to-next-to- leading-order calculation of the charm contribution to the charged decay ) 
and for rare leptonic B decays. (The comment "meas. (or lattice value) not in fit" in some rows 
indicates that for the corresponding determination the direct measurement of the quantity (or 
the LQCD prediction) has been excluded from the fit.) 

Finally, Table ^ gives the standard CKM fit results for the absolute values of the CKM 
elements and the magnitudes and phases of the sides of the non-rescaled UT. The results of the 
fit are shown in the (p, 7y) plane in Fig. The outer contour of the combined fit corresponds to 
95% C.L. exclusion. Also shown are the 95% C.L. regions for the individual constraints entering 
the fit (the constraint from i?+ t^v^- is not shown, although it is included in the fit). 

The global CKM fit performed so far contains all relevant information collected by the ex- 




Figure 8: Confidence levels in the {p, rj) plane for the global CKM fit. The shaded areas indicate 



95% C.L. allowed regions ^51]. 





Figure 9: Confidence levels in the {p,r]) plane for global CKM fits using only tree-level (left) and 
loop-induced (right) inputs. The shaded areas indicate 95% C.L. allowed regions. 



periments. From the new physics perspective, it is interesting to confront the measurements 
according to their sensitivity to new physics contributions. Figure ^ shows on the left plot the 
constraints that originate from mainly tree-level processes, together with their combined fit. 
The right plot shows the constraints from loop-induced processes. To fix the length scale of the 
UT and the constraints on A and A from the tree-level determinations of the CKM elements 
{Vudl, {Vusl and \ Vcb\ are used. If 7 is extracted from the measurement of a using /3 from mixing- 
induced CP violation as input, it is effectively a tree- level quantity, because the isospin analysis 
isolates the AI = 3/2 component in the decay amplitude, which is assumed to be standard (this 
corresponds to the SM4FC case defined in Ref. |13|). Consequently, the constraint for 7 that 
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enters the tree-level plot is the average of the direct measurement of 7 via open-charm processes, 
and the value obtained from ir — a — f3, from which new physics in mixing cancels. This pro- 
vides the first determination of (p,??) from (effectively) tree-level processes. Good agreement is 
observed between the tree-level and loop-induced constrained fits. 



7 The CKM Matrix and the Search for New Physics 

7.1 New Physics Extensions of the Standard Model and their Relation to 
CKM 

If the constraints of the SM are lifted, K, B, and D decays are described by many more 
parameters than just the four CKM parameters and the W, Z and quark masses. The most 
general effective Lagrangian at lowest order contains approximately one hundred flavor-changing 
operators, and the observable effects of any short-distance, high-energy physics are encoded 
in their Wilson coefficients. Therefore, all measurements sensitive to different short-distance 
physics are important, and overconstraining the magnitudes and phases of CKM elements could 
detect deviations from the SM. For example, Am^, T{B — > pj), and T{B Xd£~^i~) are all 
proportional to iVt^V^^p in the SM, however, they may receive unrelated contributions from new 
physics. Similar to the measurements of sin 2/3 from tree- and loop-dominated modes, testing 
such correlations, and not simply the unitarity of the CKM matrix, provides the best sensitivity 
to new physics. 

In the kaon sector both CP-violating observables, £k and e'j^, are tiny, so models in which 
all sources of CP violation are small are viable. Owing to the measurement of sin2/3 we know 
that CP violation is an 0(1) effect and only the flavor mixing is suppressed between the three 
quark generations. Hence, many models with spontaneous CP violation that predicted small 
CP-violating effects are excluded. However, model-independent statements for the constraints 
imposed by the CKM measurements on new physics are hard to make because most models that 
allow new flavor physics contain a large number of new parameters. 

There are large classes of models in which the biggest effects of new physics occur in transitions 
involving second- and third-generation fermions, thus escaping the bounds from kaon physics 
(see, e.g., Refs. [|165|| and |166(| ). For example, in SUSY GUTs, the observed near-maximal 
i^fj, — '^T mixing may imply large mixing between sn and b^, and between their super symmetric 
partners sji and bji [167]. Although the mixing of the right-handed quarks does not enter the 



CKM matrix, the mixing among the right-handed squarks is physical. It may entail sizable 
deviations of the time-dependent CP asymmetries in penguin-dominated modes from sin2/3 
discussed in the next section, and an enhancement, compared with the SM value, of Am^. 
Another recently popular scenario with large effects in flavor physics involves Higgs-mediated 
FCNCs in the large tan /? region of SUSY (tan (3 is the ratio of the vacuum expectation values of 
the two Higgs doublets, and tan/3 ~ mt/rrih is motivated by the unification of top and bottom 
Yukawa couplings predicted by some GUT models). Neutral Higgs-mediated FCNC couplings 
that occur at one loop increase rapidly with tan (3 and can exceed other SUSY contributions by 
more than an order of magnitude, giving an enhancement of B{Bd^s ~^ A*^/" ) proportional to 
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[168|. This scenario also predicts a suppression of Anig proportional to tan^/3 |169|. As 
discussed below (see Fig. B^B^ mixing may still receive a significant positive or negative 
new physics contribution even after the measurement of Anig. 

The existing data on rare B decays set severe constraints on (SUSY) model building. In 
particular B{B ^ ^s7) provides strong parameter bounds. Both the SUSY and the SM con- 
tributions enter at the one-loop level, and the charged Higgs diagrams always enhance the rate, 
whereas the chargino contributions can reduce it if /i > (/u is the coefficient of the HuHd 
term in the superpotential, and Hu^d are the Higgs fields that couple to the up- and down-type 
fermions). One reason that B — > Xg^ is the only B decay that enters the bounds in many 
analyses (e.g., the constraints on dark matter candidates |17(]| ], usually the lightest neutralino) 
is because such studies often assume the constrained MSSM (MSUGRA), which implies minimal 
flavor violation [171], and relate the charged Higgs and the chargino contributions to each other. 
The B — > X<j7 rate can be enhanced significantly even without new sources of fiavor violation, 
i.e., with all flavor changing interactions proportional to CKM matrix elements. Nontrivial fla- 
vor physics due to TeV-scale new physics would indicate deviations from this parameterization, 
which could weaken the bounds from B — > Xg^. Then the constraints can become more complex, 
and other B decays would also be important in constraining the new physics parameters. 

As mentioned above, not all CP-violating measurements can be interpreted as constraints on 
the (p, ??) plane. For example, probing with B^ J/ip4' whether Ps, the order angle of the 
squashed UT defined in Eq. (21), agrees with its SM prediction is an important independent test 
of the theory. Other comparisons between the SM correlations in B and K physics can come 
from future measurements of — > tt^uT' and — > tt~^u17. These loop-induced decays are 
sensitive to new physics and will allow a determination of /? independent of its value measured 
in B decays [ll72| . 



7.2 CP Asymmetries in Loop-Dominated Modes 

The FCNC 6 — > s transition is mediated by penguin diagrams. It can have any up-type quark 
in the loop, so its amplitude can be written as 

Ab^s = mtVtbVts + mcVcbV*s + muVubV*^ 

= (m, - mt)VcbV:g + {niu - mj)^^!/*, = 0{X^) + O(A^) , (69) 

where the unitarity of the CKM matrix has been used in the second step. In the SM, the ampli- 
tude is dominated by the first, VcbV*g, term, which has the same weak phase as the amplitude in 
B^ — > J/ipK^ decay. We expect HA/AI — 1| = 0{X'^), and the time-dependent CP asymmetry 
parameters are given to a similar accuracy by Sb^sqq ~ —r]f^psin2P and Cb^sqq ~ 0. 

Owing to the large mass scale of the virtual particles that can occur in the loops, additional 
diagrams from physics beyond the SM, with heavy particles in the loops, may contribute. The 
measurement of CP violation in these channels and the comparison with the P-to-charmonium 
reference value is therefore a sensitive probe for physics beyond the SM. A discrepancy between 
Sb-^sqq and sin2/3 can provide an indication of new physics. If the SM and new physics contri- 
butions are both significant, the deviations of the CP asymmetries from sin2/? become mode 
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dependent, because they depend on the relative size and phase of the contributing amphtudes, 
which are determined by the quantum numbers of the new physics and by strong interactions. 

The important question is how weh can we bound the contribution of the second, CKM- 
suppressed, term to the b — > sqq transition in Eq. (p9|)? This term has a different weak phase 
than the dominant first term, so its impact on Sb^sqq depends on both its magnitude and relative 
strong phase. Naive factorization suggests that for q = s the suppression of the second term 
is likely to hold because it would require an enhancement of rescattering effects to upset this. 
However, for q = u, there is a color-suppressed h ^ u tree diagram, which has a different weak 
(and possibly strong) phase than the leading penguin amplitude. For q = d, any light neutral 
meson formed from dd also has a uu component, and there is "tree pollution" again. The 
decays to tt^K^ and ujK^ belong to this category. The mesons rj' and /o(980) have significant ss 
components, which may reduce the tree pollution. Neglecting rescattering, the three-body final 
state K°K^K° (reconstructed as K'^,K^K'^) has no tree pollution, whereas — > K^K 
(excluding (l)K^) does. 

Recently QCD factorization Il7|,|l7| and SCET |l75| was used to calculate the deviations 
from sin2/3 in some of the two-body penguin modes. It was found that the deviations are the 
smallest (< 0.02) for 4>K^ and r]' . This is fortunate because these are also the modes in 
which the experimental errors are the smallest. The SM shifts enhance —rjfSf (except for pK^) 



using 173 , 174|, while suppress Sn'Ks using jl75| ]. SU(3) flavor symmetry has also been used to 
bound the SM-induced deviations from sin2/3 [ |176| , |177|| . Owing to the lack of information on 
strong phases and the weak experimental bounds on some b — > dqq mediated rates, the resulting 
bounds tend to be weak. An exception is n^Kg, where SU(3) relates the relevant amplitudes 
to Tr'^Tr'^ and K~^K~ |177]. The theoretical understanding of factorization in three-body decays 
does not yet allow accurate bounds on —rjfSf — sin2/3 to be computed. 

There has been considerable excitement about these measurements in the past few years, 
which has been somewhat damped by recent updates from the experiments. The results for the 
measurement of sin2/3efr (which is equal to —ijfSf for the final states with unique CP content) 
from the various penguin modes are compiled in Fig. 10. If one restricts the modes to those with 
the potentially smallest theoretical uncertainties, i.e., the final states (l)K^, r]'K^, and K^K^K^, 
and attempts to average the effective sin2/3 results, (sin2/?cfT) = 0.50 it 0.08, which is within 
2.2(T reach of the charmonium reference value. Beginning of 2005, this s-penguin average was 
0.40 lb 0.09, and because the charmonium result was larger at that time, the discrepancy between 
the sin2/3 numbers was at the 3.2o" level, which explains the popularity of the results. Better 
statistics are required to clarify the situation.^'' 

Another interesting measurement in the penguin sector is the time-dependent CP asymmetry 

in B ^ K^Tr^{K*^)^. BABAR and Belle have measured it exclusively and inclusively, with the 

averages SK*-y = -0.13 ± 0.32 and Sk^^o^ = 0.00 ± 0.28 ll79|,||. Although the B Xsj rate 

is correctly predicted by the SM at the 10% level, that measurement sums over the rates to left- 

and right-handed photons, and their ratio is also sensitive to new physics. In the SM, 6-quarks 

^^It is an interesting outcome of the QCD factorization calculations that the SM deviations from sin2/3 should 
be positive in all modes, except for p^Tf" [17J , 174 1. Taking the expected deviations into account would enhance 
the significance of the current discrepancy between the charmonium-based and the penguin-based results. 
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Figure 10: Comparison of world average s'm2(3(^cS) results from penguin- dominated decays, and 
the charmonium reference value [178, 12^ \12^ , [^^. The light shaded band shows the average 
of the results from (t^K^ , rj'K^, K^K^Kg, considered to be theoretically cleaner than the 

other modes. The result for ir^ir^Kg is not shown, because it has a much larger error. 



mainly decay to S'Jl and 6-quarks to S7i?, so their interference is suppressed, proportional to 
rj^ = A{B^ — > X f^jji)/A{B^ Xj^^l). If only the electromagnetic penguin operator, 0-j ~ 
's a'^'^ F^y{mi,PR + msPL)b contributed to the rate, it would give Sx'-y = —'^{'ms/mi,) sin2/3 |180|| . 
This also holds in the nonresonant B — > K^tt^^ case |181|. Grinstein et al. ]182[ recently realized 
that four-quark operators contribute to r that are not suppressed by ms/m^,. The numerically 
dominant term is due to the matrix element of O2 = {c'j^PLb)(s ^^P^c), and its contribution 
to the inclusive rate can be calculated reliably, T{B^ Xs-fR)/T{B° ^sIl) ~ 0.01 p^ . 
This suggests that for most final states, on average, r ~ 0.1 should be expected. A SCET 
analysis of the exclusive decay proved the power suppression of the amplitude ratio, A{B^ — > 
K°*-fR)/A{B^ K^*-iL) = 0[{C2/2,C^){AQCD/mb)] ~ 0.1 ll8|, but the uncertainties are 
sizable. 



7.3 The D System 

The D meson system is complementary to K and B mesons, because it is the only neutral-meson 
system in which mixing and rare decays are generated by down-type quarks in the SM (or up- 
type squarks in SUSY). Therefore, flavor and CP violation are suppressed by both the GIM 
mechanism and by the Cabibbo angle. As a result, CP violation in D decays, rare D decays, 
and D^D^ mixing are predicted to be small in the SM and have not yet been observed. There 
are well- motivated new physics scenarios, such as those based on quark-squark alignment |183(| , 
that predict D^D^ mixing to be of order A^, so that if the current bounds improve by a factor 
of a few, these models become fine tuned. 
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The most interesting experimental hint for a possible D^D^ mixing signal so far is the lifetime 
difference between the CP-even and CP-odd states 11184 



r(CP-even) - r(CP-odd) ,^,,,,,0/ 
r(CP-even) -|- r(CP-odd) 

Unfortunately, this central value alone cannot not be interpreted as a sign of new physics, because 



of possible long-distance hadronic effects [185]. At the present level of sensitivity, CP violation in 
mixing or decays, or enhanced rare decays, such as D — > TTi~^i~ would be the only clean signal of 
new physics in the D sector. Because the i?-factories are also charm factories (however, without 
coherent production of D^D^ pairs), it is possible to observe such unambiguous signals of new 
physics every time the analyzed data sample increases, independent of the tight constraints in 
the K and B sectors. 

7.4 The System 

Although all present results show agreement with the SM, it is still possible that the new physics 
flavor structure is quite different from that in the SM, with 0(1) effects in sectors of flavor physics 
that have not yet been probed with good precision. As mentioned above, many well-motivated 
models can lead to large new physics effects between the second and the third generations, while 
leaving the flavor-changing transitions between the other families unaffected. The sector, 
which is currently being investigated at the Tevatron (cf. Sec. |3^) and will soon be studied by 
the LHC experiments, is suited to test such scenarios. 

Compared with the B^ mesons, the properties of the B^ are still less well known, although 
the measurement of Am^ p5|,|96t strongly indicates that B^B^ mixing is SM-like. Although the 
uncertainty of Am^ is already at the few percent level, we have no significant information yet on 
the BgB^ mixing phase, /3s, and the width difference, AF^, is poorly known. Many extensions 
of the SM predicted large modifications of Am^ owing to contributions from new particles to 
B^ oscillations. Such models, or such parameter regions of models are now ruled out. The LHC 
experiments will determine sin2/3s through the time-dependent measurement of mixing-induced 
CP violation in the decay B^ J/i^'P (and similar final states). The expected precision 
of better than 0.02 for the combined LHC experiments will barely be enough for a significant 
measurement if S^p^j/^^ is as predicted by the SM. However, for the purpose of the search for 
new physics, which could still occur at the C(l) level, only the absolute error on sm2f3s matters. 

^^During the Summer of 2005, the BeUe experiment performed a three-day engineering run comprising a five- 
point energy scan between 10.825 GeV and 10.905 GeV, to locate the T{5S) peak, and data taking at the T{5S) 
(10.869 GeV). A data sample corresponding to an integrated luminosity of 1.9 fb~^ was collected during that 
period ll86|. T{5S) resonant events decay into pairs of neutral and charged i?'-*-' mesons as well as neutral B^*^ 



mesons with opposite flavors. Belle found a fraction of events with bI*' pair production over the full bb production 
(including the continuum) at the T{5S) of (16.4 ± 1.4 ± 4.1)%. Among these, T{5S) B'^Bl events dominate. 
All inclusive and exclusive measurements performed are found in agreement with earlier results on the T{5S), 



with less statistics, obtained at the CLEO experiment 187 1. 

^^As all pseudoscalar decays to a pair of vector mesons, such as ^ J/tpK*'^ or _B" , the final state 

in Bs — > J/ip4> is an admixture of CP-even and CP-odd components. To avoid a dilution of the CP-asymmetry 
measurement, they must be separated by an angular analysis, which complicates the measurement. 
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Figure 11: Confidence levels of the hd,crd (left) and hs^Ug (right) new physics parameters. The 
dark, medium, and light shaded areas indicate 10%, 68% and 95% C.L. allowed regions, respec- 
tively. The SM corresponds to hd = hg = 0. 

Because this phase is suppressed, it can be predicted accurately (in terms of its absolute 
error) by the global CKM fit, giving sin2/3s = 0.0365 ± 0.0021 [ll8|,|ll. Hence this represents 
the first precision measurement, sensitive to a new physics phase, that can be confronted with a 
precision prediction. It could become the most sensitive test of the KM theory for CP violation. 

The class of rare decays that proceed through penguin diagrams, such as i? — > K*^ and 
B K*i~^i~, will be augmented by the system through (among others) the measurement 
of the decay B^ A*^/^^ at the LHC. The best current limit on the branching fraction is 
8 X 10"^ at 90% C.L. |l8|], which is well above the SM model expectation of 3.4 x 10"^ (whose 
dominant theoretical uncertainty is due to f'^J |19C(| . Expectations for the LHC experiments 
are hard to estimate owing to the difficulty in fully simulating sufficient background statistics, 
but all experiments expect to see a SM signal within one year. The decay B^ A^^A*" will be a 
further challenge due to its even lower branching fraction in the SM. It should be observed by 
all experiments after several years of data taking. A precision measurement of B{B^ ^'^ ^~)/ 
B{B^ /i+/x~), which has a clean relation to Arus/ Am^, will be beyond the scope of the LHC. 

7.5 Model-Independent Constraints on New Physics in B Oscillation 

In a large class of models the dominant new physics effects are new contributions to the B^B^ 
and BgBg mixing amplitudes |191[|, which can be parameterized as 



M 



12 



Mf2^(l + /i,,,e^-''-) 



(71) 



The allowed regions for ad vs. hd and ag vs. hs are shown in Fig. 11 (left and right plots 



respectively). Numerically, one finds hd = 0.23 _o'23; corresponding to the 95% C.L. upper limit 
hd < 1.7 [|3|. The constraints on the new physics parameters are reasonably tight because (p,rj) 
has been determined from (effectively) tree-level measurements independent of new physics in 
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B^B^ mixing (before the a and 7 measurements most of the plotted region was allowed) 
192,193,194]. In B^B^ mixing, the Im.e^^'^'i > region is more strongly constrained than the 



Im €"^'"1 < region due to the mild tension between sin 2/3 and \Vub\- Although the constraints 
are significant, new physics with a generic weak phase may still contribute to the B^B^ mixing 



amplitude of the order of 30% of the SM. Because Eq. (71) has two new parameters that need 
to be determined from the data, measurements irrelevant for the SM CKM fit, such as the CP 
asymmetries in semileptonic B^ and B^ decays become important for these constraints [ 188 , 195 , 



19^ ; the former is measured as^o A^^ = -(0.5 ± 5.5) x lO'^ ^T^. In B^B^ mixing, even after 



the measurement of Anis, new physics and the SM may still give comparable contributions to 
the mixing amplitude, since we have no information on the mixing phase yet [166, 195 [. This 
will change as soon as there is data on the CP asymmetry in B^ J/ip4''i even a measurement 
with an uncertainty as large as 6(SBs^i)<l>) ~ ^-^ would bring the constraints on hs,as to a level 
comparable to those on h^,ad |195([ . At present, the lifetime difference between the CP-even 
and odd B^ states gives a useful additional constraint [197, 19^ , 195|[ (not included in Fig. \n\ ), 



however, this is largely because the central value disfavors any deviation from the SM. After 
a nominal year of LHC data this constraint will probably be less important due to theoretical 



uncertainties. Similar results for the constraints in K^K^ mixing can be found in Refs. [166, 19 



8 Conclusions and Perspectives 

The i?-factories and the Tevatron experiments have provided a spectacular quantitative confir- 
mation of the three-generation CKM picture. Even more interesting than the improved determi- 
nation of the CKM elements is the emergence of redundant measurements that overconstrain the 
CKM matrix. The measurements of CP asymmetries, mixing, and semileptonic and rare decays 
allow us to severely constrict the magnitudes and phases of possible new physics contributions 
to flavor-changing interactions in a variety of neutral- and charged-current processes. We can 
meaningfully constrain flavor models with more parameters than the SM; in particular, within 
the last two years, the comparison between tree- and loop-dominated measurements significantly 
bounds new physics in BB mixing. This illustrates that it is the multitude of overconstraining 
measurements and their correlations (and not any single measurement) that carry the most 



compelling information. Nevertheless, as can be seen in Fig. 11, new physics in B^B^ and in 
B^Bg mixing may still be comparable to the SM contribution. 

Having analyzed these results, one may ask where flavor physics should go from here. The 
sensitivity to short-distance physics will not be limited by hadronic physics in many interesting 
processes for a long time to come. The existing measurements could have shown deviations from 
the SM, and if there are new particles at the TeV scale, new flavor physics could show up any 
time as the measurements improve. If new physics is seen in flavor physics, we will want to study 
it in as many different processes as possible. If new physics is not seen in flavor physics, it is still 
interesting to achieve the experimentally and theoretically best possible limits. If new particles 
are observed at the LHC, the measurements from the flavor sector, as well as flavor-diagonal 
To avoid contamination from iSf mesons, only the measurements performed at the T(45') are included in this 



20 



average [197|. In the SM ylsL ^sli so one may average all data to bound \q/p\ — 1, quoted on p. ^ 
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constraints from electric dipole moments, will strongly confine the degrees of freedom of the 
underlying new physics. 

The large number of impressive new results speak for themselves, so that it is easy to sum- 
marize the main lessons we have learned so far: 

— sin2/3 = 0.687 it 0.032 implies that the B and K constraints are consistent, and the 
KM phase is the dominant source of CP violation in flavor-changing processes at the 
electroweak scale. 

— The difference between the measurements of sin2/? from B decays to charmonium plus 
and from the penguin-dominated rj'K^ mode (or the average of (f)K^, rj'K^ and K^K^Kg) is 
approximately 2a at present. More data are required to make this comparison conclusive. 

— The measurements of a = (lOOlg^)" and 7 = (62125)° are in agreement with the SM 
expectations, and begin to improve the constraint on p, r/. They also provide significant 
bounds on new physics in B^B^ mixing. 

— The CDF measurement of Am^ = 17.31 1q'^| ib 0.07 ps^^ (in agreement with the SM 
expectation, 21.7^4 2) together with the results of the other experiments lift the evidence 
for BgBg oscillation to the 4cr level. 

— The determinations of \Vcb\ = (41.5±0.7) x 10"^, mb{mb) = 4.18±0.04 GeV, and mc(mc) = 
1.22 ±0.06 GeV from inclusive semileptonic B decays reached unprecedented precision and 
robustness, as all hadronic inputs were simultaneously determined from data. 

— Strong experimental evidence for the leptonic decay t'^u-j- has been found, with a 
world-average branching fraction of (10.4 12;?) x 10~^, in excellent agreement with the SM. 

— Ax+^- = —0.115 lb 0.018 implies large direct CP violation disproving "B-superweak" 
models. It also signifies that sizable strong phases occur in this decay. Many more decays 
exhibiting direct CP violation will emerge from the increasing data samples. 

— The list is much longer, e.g., improvements in the determination of \Vub\', observation of 
B —>■ Xsi'^£~ and B — > decays, and a zoology of new Dg and cc states. 

The next few years promise similarly interesting results (in arbitrary order): 

— One of the most important outstanding questions is the comparison of sin2/3 from tree 
and penguin decays. Three times more data than currently used in the analyses could be 
accumulated until data taking of the first generation i?-factories ends. 

— A measurement of sin2/3s with mixing-induced CP-violation analysis of B^ J/4'4' (and 
similar) decays. If sin2/3<j can be measured at the 0.02 level (absolute error) at the LHC, 
it will constitute the most sensitive probe so far for new CP-violating phases. Together 
with other SM-like observations in the K and B sectors, it would suggest that the flavor 
sector of any new physics at the TeV scale approximates minimal flavor violation. 

— Improvements in the determination of a and 7 can be expected. The angle 7 will also be 
measured by LHCb. This will further constrain the right side of the UT, serving as the 
SM reference for the ratio Arud/ Amg. Reduction of LQCD errors is also necessary for this 
comparison. 
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— The rigor of the determination of |V^(,| may approach that of \ Vcb\, and its error wih also be 
reduced. It determines the side of the UT opposite to f3, so any progress directly improves 
the accuracy of CKM tests (the error with continuum methods reach an asymptote at 
approximately 5%). 

— The experimental precision approaching the limits from the theoretical uncertainties in 
B Xgj and B Xst^i~ (for branching fractions, CP-violating asymmetries and 
forward- backward asymmetries) will impact model building. The (continuum) theory is 
most precise for inclusive decays, which, however, cannot be measured at the LHC. 

— Improving the precision of the current Sk*^ = — 0.13ib0.32 measurement will be important 
to constrain certain extensions of the SM. 

— Improving the precision of the current Asl = — 2(|g/p| — 1) = —(0.5 ib 5.5) x 10~^ mea- 
surement is also important, in addition to a measurement of in B^B^ oscillation. 

— Future studies should firmly establish and precisely measure the B ^ and B ^ tv 
rates, and try to measure B X^'J inclusively. They should also search for B^^g — > i^i^ 
and other rare or forbidden decays. 

— Similar to the "new" cs and cc states discovered by the S-factories, new physics could also 
be discovered in the charm sector. With increasing sensitivity, one may find clear signs of 
new physics in D ^ irt^l" or through CP violation in D^D^ mixing. 
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